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Abstract—The momentum transfer in a two-phase stream, consisting of particles carried by a 
fluid, was studied by applying the statistical theory of turbulence. Statistical relations showed 
that the correlation, scale and intensity of one phase can be calculated from those of the other 
phase. Approximate solution indicates the basic parameters affecting momentum transfer. 


Résumé—En appliquant la théorie statistique de la turbulence, l’auteur étudie le transfert 
des quantités de mouvement dans un courant 4 deux phases constitué de particules transportées 
par un fluide. 

Des relations statisques montrent que les caractéristiques, proportion et intensité d'une 
phase, peuvent se calculer a partir des caracteristiques de l'autre phase 

A l'aide de solutions approchées, l'auteur trouve les paramétres de base, afféctant le transfert 


de la quantité de mouvement 


INTRODUCTION conveying, efforts have been made toward finding 
Tue characteristics of momentum transfer in a the relation between the motion of solid particles 


; “<a i" - - 
two-phase stream consisting of liquid or solid and the fluid stream (4, 5). 


particles in a gas, or solid particles in a liquid, is T) 


purpose of this pape is to present the 


of fundamental interest in the following fields - ; : 
relations between the intensities and scales of the 


of study and techniques of measurement , 
. two phase S arrived at from the statistical theory 
} eo 4 ; . : 

1. Inusing solid particles astracersinstudying of turbulence, together 1 an approximate 
flow of fluids, the proximity of the particle motion gojytion which indicates the basi parameters 


ii 


> » bk | ** ' 
to stream motion should be known [1]. affecting the momentum transfer between the 


2. Ina fluidized bed of catalyst, the diffusivity two phases 
of the continuous phase can be measured with The two phase fluid under consideration will 
relative ease, but the diffusivity of the particles be a fluid stream carrying spherical solid particles 
has not been taken care of [2]. of uniform size. In this system, mass transfer is 
absent and heat transfer is negligible between 


8. In the combustion of solid fuel particles or 


liquid fuel droplets, the diffusion between the the two phases. The following assumptions have 
fuel and air is influenced by the relative motion been made in the following analysis : 


of the fuel particles and the air stream [3}. The turbulence of the fluid stream is isotropic 


4. In sedimentation studies and pneumatic and non-decaying. 
* Research project sponsored by the Research Corporation, New York, N.Y. 
+ Assistant Professor of Mechnical Engineering, Princeton University, Princeton, N.J. 


** Numbers in brackets refer to references at the end of the paper. 
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The relative velocity between the solid parti- 
cles and the stream gives, at all times, a 
Reynolds number of the particies of order 
1 or less. 


The particles are relatively far apart so that 


the inter-action between the particles can b 
neglected (see Appendix 1). 

The average velocity of the solid particles is 
equal to that of the fluid ind 


constant.* 


stre is a 


The case of pendulum motion of a | 
turbulent stream was studied by [6] 
However, application of the Eulerian rrelation 
in the present case would require the ision of 
the probability of the existence of a | 
given [4); 
introduced here so that the velocity 


location Lagrangian correlation is 


ia single 
In the 
ition of 


té¢ red by 


particle is described along its jour! 
Lagrangian case, the rapidity of fluctu 
veloc ity in time or its character Is not 
the presence of the average velocity the main 
itistical 


differ 


the 


stream. Therefore, the Lagranyian s! 
frequency function f(y) is fundament 
that of the 


difference is pointed out in the follow 


ent from Eulerian, { 
the 


ated to thi 


In adopting th Lagrangian cor ition, 


velocity of a fluid particle becomes r 
or its velocity f its 


time co-ordinate wing 
motion, and the velocity of a fluid particle at a 
the 
time 
the 


re spect to 


fixed point becomes insignificant ugh 


tive 


relation between the space location a 
taken of thr 


integral of velocity components with 


coordinates can be care ugh 


ym in the 


and the 


time). Following the fourth assumpts 


above, the velocity of the solid partick 
velocity of the stream are related t 


thea 
it 


Sanne 


SA Tre 


time co-ordinate, although not the space 


In other words, for thy under 
consideration, the of a 
particle and the velocity of fluid motion are 
ind 


ist 


co-ordinate. 
solid 


velocity single 


related to the same time coordinate, we are 
concerned only with the velocity of the fluid 


stream in the neighbourhood of a solid particle. 


* The problem of acceleration of the particles from rest 
by a fluid stream of uniform velocity was sufficiently 


treated in reference [5). 


58 


FUNDAMENTAL RELATIONS 


Here the concept of spectrum of t irbulence 
of 
stream, which flows horziontally 
UU, the 

f the st 


ol turbule nee ol 


made us¢ Due to isotropy of turbulen 


with a 
vellocity iXial 
to l 


representative of all three compo 
of the 


average omponent, 


(relative 


velocity spectrum 


A 


and n is the Irequency. 


is the amplitude or th 


The 


whe re = 


intensity 


ke nce ; 


Lim , 


where f.(n) is the Lagrangian statistical frequency 


function, and 


. , 
i tn) dn 


KE xte nding Tay LOR’S hyp 


system : 


=» & 


‘ 


. u(t) ull 
R. 
u* 


which is a time correlation rather than that of 


equating Eulerian space correlation to Eulerian 


time correlation in the Eulerian system of 


turbulence spectrum. 
By the definition of f,(n), R, can be expressed 
in terms of a Fourier integral : 


ny 


Vu? 


f(n)cos2an Edn { (n)cos2a 


+ 


v0 0 


where dy / u?® dl, and for non-decaying tur- 


Vv u* 


Rand fi(n) are 


bulence, » Ilence, as analogous to the 


Eulerian cacs, Fourier trans- 
forms [8]. 

The above relation implies that the velocity of 
a single particle of the stream is described along 
its journey, and the rapidity of its fluctuation in 


time is not altered by the average velocity U of 
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the main stream. The concept of diffusion by 
continuous movement gives [9] : 


(X?) = a R, dé 


where A® is the average displacement. 
Further, the scale of turbulence of the stream 
is given by : 


T 


l, = | 8a, 


0 


(6) 


and the A, defined by the parabola in Fig. 1 is 
given as: 


(7) 
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Fic. 1. Approximations to Lagrangian correlation curve. 


The case of motion of a spherical pendulum in 
a viscous fluid was treated by Stokes as pre- 
sented in Reference [10]. Since in this present 
study we are interested in small solid particles 
having very small motion relative to the stream 
(in terms of relative frequency between the parti- 
cle and the stream), the correction for increase in 
apparent inertia and the diminution of oscillation 
[10] can be neglected. Therefore, neglecting the 
squares and products of the relative velocity 
between the fluid and the solid particle, and for 
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constant average velocity of the main stream,* 
the motion of the particles can be described by 
(for small relative motion defined by assumption 
(2) in the above) : 


d Uy 


(8) 


18 


F 
d,* Pp 


[12] 

F has the dimension of (time) and is character- 
istic of the time response of the particle. Equation 
(8) was obtained through simplification by the 
comparison of order of magnitudes rather than 


simply an engineering postulate. Solving Equa- 


| 


(9) 
2mn 
F) 
To calculate the vertical component of the tur- 
bulent velocity of the particle, neglecting the 
gravity effect on the stream but considering the 


tion (8) with u given in Equation (1) : 


a] 


a 
— sin 2mrnt — cos 2mrnt 
2mn 


(2nt — 6,), 


where @,, is the phase lag and 6, = tan™ 


gravity effect on the particle : 


F(v —v,) +4, (10) 


where g is’ the acceleration due to gravity, which 
gives : 


2nn A 
— ante: | 
0p = FL aT am) 

E sin 27nt — cos 2ent| - z (11) 


The intensity of turbulence of the particles are 
respectively : 


(12) 


* The case of an accelerated field without turbulence 
was treated recently by Grteert and others [11]. 





(18) 


Therefore, effect is 
considerable, the turbulent motion of the particles 
For excessively large particles such 


in general, where gravity 
is anisotropic. 
that the terminal velocity is much greater than 
u2, vp’ has no meaning, nor u,*, as the particle 
would fall quickly to the bottom of the duct. (See 
Appendix “or carrying capacity of a turbulent 
stream.) For cases of small particles and there- 
isotropy of the turbulent motion 


Hence 


fore large F, 
of the particles can be assumed. 


2nn 

F) 
By the definition of the frequency 
(Equation 2) : 


: (14) 


function 


l 
A? =v dn, 
@ Ae u* f(n) dn 


wf, . 


9 
1+ (7) ' 


Equation (15) can 
However, f(n), the 


If function f,(n) is known, 
be used to calculate u,’. 
Lagrangian frequency function, is not as well- 
known as that of the Eulerian [18]. In the 
following, an approximation was made from 
R,, which is measurable. 


The relative velocity between the particle and 


and the stream can be represented by : 
Au = u — u,. (16) 


The intensity of turbulence of the particles 
relative to the stream can be obtained from : 


Introducing Eqs. (1), (9) and (16) into Eq. (17), 
it can be shown that, without any further simplifi- 
cation, 

T 


(Au)*® dt = 


Lim 
T—+@® 2T , 
T 


4° 


The Lagrangian correlation of the particle 
motion can be shown to be : 


R, — “elt) p(t €) _ 





In terms of y, substituting »/+/w? for ¢ in Eq. (19) 
and postulating as in Eq. (15), the Lagrangian 
correlation of the particle can be expressed in 
terms of the Fourier integral : 


— fA) cos Ben Fann 
1+ (FF) 





Ri() = (20) 


Sln) 
' 2an\2 
+ (F) 
0 
The scale of turbulence of the particles : 
T 


hi = { Ry ay, 


Ry (1), 


dn 


R, (n) = 
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where 


so B 


Therefore, knowing the Lagrangian correlation 
(scale) and intensity of one component of the 
two phases, and the properties of the stream and 
particle (F’), the scale and intensity of turbulence 
and the correlation of the other phase can be 
calculated. 

Other relations of interest are : 

(a) The viscous dissipation due to each particle 
of mass m can be calculated from: (for three 
components) 


n) dy * j Ri(n) dy’. (22) 


T 
l 
= 8 Lim oF 


T+@ 


Ww 


» (drag force) (Aw) dt 


z 
(Au*) dé, 
-f 


=8 Fm Lim 35 


Drag force + Fm(Au). From Eq. (18), 
W, = Om d, (u® — u,*) (28) 


(6) The ratio of particle diffusivity to eddy 
diffusivity, from Einstein's equation of diffusivity, 
is [2) : 


for : 


(24) 


The evaluation of the integrals obtained in the 
above from the curve of f,() analogous to /(n) 
given in Reference [18] is tedious. The relative 
magnitude of various quantities and basic 
parameters affecting momentum transfer can be 
seen from the following approximate solution : 


APPROXIMATE SOLUTION 


The curve of 
R= ol" 


has the same curvature at the apex as the parabola 


(25) 


7 
R,=1-— (26) 


A,* 


where A, was explained by TayLor as a measure 


of smallest eddies responsible for the dissipation 
of energy. This parabola, in many cases, is a 
sufficient approximation of R,. 

The curve of R,(y) from experiment by TayLor 
is reproduced in Fig. 1 together with curves given 
in Eqs. (25), (26) for comparison. Use of corre- 
lation curve obtained in Reference [2] showed 
similar type of agreement. 

Using the approximation given by Eq. (25), 
from Eq. (5) : 


2 a/m A, Pie liad bal 


Vu 


f,{”) , (27) 


From eq. (15), 


_ 2 
us? = 


7 any 


= lie 
a K+ S 3 


ere.. i 
9¢ 


or 


u,* =u GK); G(K) <1; (29) 


* It would be interesting to notice that, with the 
Eulerian system, using the approximation : 
R, ~_ e/a 
The corresponding frequency function is : 


fin) = 224 parermier, 


7 Ye A. (b) 


The difference between /,(n) and f/(m) can be seen by 
comparing eqs. (27) and (a). Further, TayLon gave : 
2 A® 
-=—-— « ] 


9 A, 
or A, = 0-5 and | 0-51, from eqs. (28) and (b). 
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Pa 
\ FA. 


1 /d, 2 d,, : l 
‘ } a | (=) (°) = 
the 


transfer between the two phases can be observed. 


30) 


Re : 


from which factors affecting momentum 


The Reynolds number Re is the Reynolds number 


held at a 


\ elocity 


of the particle if the solid partick 
uniform velocity equal to the averag 
of the Physically, this 


squared number a char 


stream. »t-mean- 


Rey nolds 


Reynolds number of the viscous drag 


IS teristic 
parted 
to the soli particle by the stream to produce 
oscillatory motion of the particle 


From eq. (20), 


1,.*a'*n'*/u* 


u* 24 nr A, 
\ uz 


¥ 
} 


COS ln 
2 


€ 
2mrn\? 
1+| 


et/4n HK, » A.) /G(K), 
where 


n (K, 2) 


n 


G(K) 


105 


t (a) (K- 


+ (3) (# 


(7) —146K* +....) 


- 


+(7 ) xe-.... Ras 


The integration formulae for the above transforms 


are present in Appendix II. From eq. (25), 


R, = R,H (K. 1) / em) (82) 


2, —nH (K, 7 


" 


+InG(K), (83) 
. 

from which A, can be solved, since u,* is known 
from [29]. However, A, can again be approxi- 
mated by finding the parabola having the same 
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curvature at the apex of the R)(») curve 


pressing R.(») in terms of : 


evn th + (2) 


bes 


and the curvature at the apex of the characteristi 


be 


l 


2 J (k) 

G (k) 

2[G(K) I (K)| 
A,? [G (k)} 


PR: 
dy 


9 


— 


parabola is : 


d* BR. 
d»* A? 
Hence 


A? 


AS 


G(K) 
I1(K)/G(kK 


terms of the above 


The 


A? 


A; 


first few expressi 


291 


16 


K* 


The scale of turbulence : 


Comparing eqs. (28), (37), and (38), J, 











- -o 


Fic. 2. Ratios of intensities, scales, and diffusivities of 


particles to those of steam at various K, 
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general. The ratio of particle diffusivity to eddy eq. (30 } 
diffusivity will simply be, from eq. (24 of the squares 


A, /u,? G (K) 
A, wu? V 1—1( K)/G(K) 


velocity between th 


3Y) becomes significant 


The values of V/ u* by*, A./A,, and E,/E fur 
various values of K are as shown in Fig. 2. Th 
relation between (Au)*/u? and K is shown in 
Fig. 3. 


Fic. 3. Ratios of intensities of relative m 
particles and steam to intensities of steam 


DIscUSSION 


The above shows that. knowing the 
scale of turbulence of one ph 
and knowing .the properties of 
particles, the intensity and s 
phase can be calculated (for sma 
between the particles and the 
Where the condition that R 
applies. measurement of onc 


whichever is more convenient mia Ti he | nenceyv oO evlecting the interact 


iny cases such as 


Functions G(A), J(K) and H(A 


ind pneuniatsi conveving can 


from integration of exact statisti mn from the experimental results in Appendix 1, 
function can be tabulated from machine calcul where average distance between the particles is 
tions for interpolation from measured quantitic greater than 10 diameters of the particle. How- 

Figs. 2 and 8 show that the above theoretical ever, in the case of fluidization, the particles 


results will not apply for large values of A might be 14 diameters apart from centre to 
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S. 


the 
interaction of the particles has to be considered. 


2 and 


centre even for 80% void (Ref. 


It is obvious that in order that the particles 
ean be suspended by a low density stream, the 
turbulence must be non-decaying and the viscous 
be 


dissipation of the particles (eq. (23)) must 


replenished by the turbulence generated from 
the boundary layer of the flow duct at the expense 
In of fluid the 
turbulence of the stream is produced by the 


i down by 


of pressure drop. case zation, 
eddies generated from the particles h 
gravity. Therefore, from the energy consideration, 
the turbulence existing in the two-pl stream 
can be calculated. An interesting problem here 
is the behaviour of the particles in the boundary 
layer of the stream. 

The the al 
applied to the following practical cases 

l. To 
pattern in the tracer technique, th 


results obtained in can be 


‘ 
. 


to flow 


get close proximity ream 
ndition is 
to choose particles of size and density so to limit 
the value of K, from estimated flow condition of 
the to the For 
instance, to get better than 98% accuracy of the 


flow pattern, K should be less than -04 as can 


stream, accuracy required. 


be seen from eq. (40), (Fig. 8). 
In 


(Au)* in eq. (28) is significant, as it represents 


9 


2. pneumatic conveying, the quantity 
the energy dissipated in suspending the solid 
particles by turbulent motion. This energy must 
be replenished by the blower besides the energy 
the 


dissipated by pipe friction 


8. 
liquid fuels, the design of the combust 


in accelerating particles and the energy 


[5] 
In the combustion of dispersed solid or 
n chamber 


It 
high 


affects greatly the diffusion of fuel particles. 
get 


~ 


; 


can be seen from the above that 
particle diffusivity, the of 
should not be too large and therefore a number 


scale turbulence 
of small combustion chambers of reasonable size 
for proper mixing should be used instead of a 
single large combustion chamber, and excessive 
atomization should not be attempted (to maintain 
> > 4 . > 
large d,/A, or large K, Fig. 2). 
4. 
although the desirability of large surface area 


In the fluidization of beads of catalyst, 


favours smal] size of beads, particle diffusivity 


L. 
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Soo 


» 
- 


an 


(Fig. 


optimum design can be achieved through the 


), however, favours large size beads ; 


use of beads of suitable size. In other words, 


large relative motion between the beads and the 
fluid stream is desirable (Fig. 3). 


CONCLUSIONS 


In a two-phase stream the characteristics of 
turbulence of one phase can be determined 


from that of the other phase. 


The fundamental parameters affecting the 


statistical properties of momentum transfer 
between the two phases are : 


d ~~ a d \ 
(=X*): (ez) (2) 


When gravity effect of the particles is signi- 


) 
rs 


p 


ficant, the turbulent motion of the particles 
is anisotropic, even though the stream is 
isotropic. The parameters affecting moment- 
um transfer will then be : 

4/ ut 


| ), (vue), (*°), (2) 

Vd, 8 eh / \GS \p, 

In general, the scale of turbulence of the 
particles is greater than that of the stream; 
the intensity of the particles is less than that 
of the and the diffusivity of the 
particles is greater than the eddy diffusivity 
of the stream, but tends to the eddy diffusivity 
of the stream. 


stream ; 


Acknowledgement—The author wishes to thank 
Dr. Ricnarp H. Hucues, Shell Development 
Company, Emeryville, California, for his many 
valuable comments and suggestions. 


APPENDIX I 
Carrying Capacity of an Isotropic Turbulent Fluid Stream— 
The approximate relation for calculating the carrying 
capacity of a turbulent stream was developed as follows : 
TAYLOR gave the energy of turbulence of a fluid stream 
as: (8) 


5 (3 wu?) 


In suspending the particles in a two-phase stream, only 
1/6 of 3 wu? is effective, when the intersection of the par- 


ticles is neglected. For a duct of cross-sectional area A 
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and length L, the energy available for suspending the 
particles will be, approximately 


pLA 
4 


wu 


Since only the components of u greater than the terminal 
velocity, V,, of the particles will be effective and 


2 
V _ d,, & (Py — P) 
: 18 


To suspend N particles, each of mass m, the energy 
required will be : 
M 
Nm V2 = (“ v2 


where M is the particle feed rate; At = 
avcrage velocity of flow. 


L 
ig and U is the 


M - p A rf 


where z is the percentage of ae greater than - f - and 
V u2 u2 

the probable value of z is given in Fig. 4, obtained from 

statistical velocity distribution developed from kinetic 

theory of gases [14). 

HaALisTrom [15] measured experimentally the maximum 
carrying capacity without build-up of a turbulent stream 
in a 8in. diameter duct at various average stream velo- 
cities, using glass beads of sizes from 105-123 microns in 
diameter. The intensity of the stream can be calculated 
from the above equation with the aid of Fig. 4: 
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Fic. 4. Probable velocity distribution. 





Measured 
Quantities Calculated Quantities 


Mean dis- 
Air flow tance 
U M rate, 


fps 


Particles between 


Vu2/U 
lbs /min | lbs /min in® particles, 
dia- 
meters 
—— 
26-2 | 0-727 

28-8 | 


1-26 
34°8 3-12 


1420 
2230 
4570 





APPENDIX II 
Integration Formula Developed for the Approximate Solution 
oO @ 
e a2n2 
- dn 
J1+0? n* , 
0 0 


Vel, Y? : 1-3-5 6° 1-3-5-7 b* 
2a 2a? ; 2 9° ; 


= G (b? /a*) - Vn 2a 


—_ ‘ 
¢ 22n2/) _ 422 4 tn * ~ 59n8 1 b9n 8... )dn 


2) 
e-azn2 ‘1 
s3°" 2 pndn = ¢-a2n2 () b?n? b*n* 
1+ b*n 

a b*®n® + b8n® - ) cos Z2pn dn 


= Fy — b* F, + o* F,- b° F, + 0° F, 
where 


Ww 
Fo = —— ¢ P*/a° 
2a 
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It is obvious that in order that the particles 
ean be suspended by a low density stream, the 
turbulence must be non decaying ana the viscous 
23 


dissipation of the particles (eq. ( must be 


replenished by the turbulence generated from 
the boundary layer of the flow duct at 
In of fluidiz 
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turbulence by 


nsideration, 


gravity. Therefore, from the energy « 
the turbulence existing in the two-pl stream 


can be calculated An interesting probiem here 


is the behaviour of the particles in the boundary 
layer of the stream 
The obtained 
applied to the following practical 
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results in the 
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can 
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the value of K, from estimated flow 
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instance, to get better than 98° 


stream, accuracy requ 


act r 
flow pattern, K should be less than -04 as 


be seen from eq. (40), (Fig. 8 
In 
(Au)? in eq. 
the 


particles by turbulent motion 


” 
-- 


pneumatic conveying quantity 
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#-—)} 


Is Significant, as 


re presents 


dissipated In suspending the solid 
This . 
be replenished by the blower besides the energy 
the 


dissipated by pipe friction | 


8. In the of 
liquid fuels, the design of the combust 


energy 


must 


nergy 


in accelerating particles and the energy 


[5 
combustion disperse 1 solid or 
chamber 


It 
high 


affects greatly the diffusion of fuel particles 
can be seen from the that 


particle diffusivity, of 


should not be too large and therefore a number 


, 


above get 


the scale turbulence 
of small combustion chambers of reasonable size 
nstead of a 


excessive 


for proper mmxing should be used 
single large combustion chamber, and 
atomization should not be attempted (to maintain 


2). 


large d,/A, or large K, Fig. 2 
4. In of beads 


although the desirability of larg 


the fluidization of catalyst, 


surface area 


favours smal] size of beads, particle diffusivity 


. Soo 
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» 


(Fig. 2), however, favours large size beads; an 
optimum design can be achieved through the 
use of beads of suitable size. In other words, 
large relative motion between the beads and the 


fluid stream is desirable (Fig. 3). 


CONCLUSIONS 


In a two-phase stream the characteristics of 
turbulence of one phase can be determined 
from that of the other phase. 

The fundamental parameters affecting the 
statistical properties of momentum transfer 


between the two phases are 


(“e*), (*), | 


When gravity effect of the particles is signi 


Pr. 


p 


ficant, the turbulent motion of the particles 
the 


isotropic. The parameters affecting moment- 


is anisotropic, even though stream is 


um transfer will then be : 
VV u* 


= ); (“2 u* *). (- P), ("* ). 
V ad,8 be lL, a 
In general, the scale of turbulence of the 
particles is greater than that of the stream; 
the intensity of the particles is less than that 
of the the 
particles is greater than the eddy diffusivity 
of the stream, but tends to the eddy diffusivity 


of the stream. 


stream; and the diffusivity of 
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ApPrpenvix I 
Carrying Capacity of an Isotropic Turbulent Fluid Stream— 
The approximate relation for calculating the carrying 
capacity of a turbulent stream was developed as follows : 
TayLor gave the energy of turbulence of a fluid stream 
as: (8) 


‘(3 u*) 


In suspending the particles in a two-phase stream, only 
1/6 of 3 wu? is effective, when the intersection of the par- 


ticles is neglected. For a duct of cross-sectional area A 





Statistical properties of momentum transfer in two-phase flow 


and length L, the energy available for suspending the 
particles will be, approximately 


Since only the components of u greater than the terminal 
velocity, V,, of the particles will be effective and 


dp 8 (Py — P) 
18 uw 


To suspend N particles, each of mass m, the energy 
required will be : 


M (At) 1,4 


i. 
Na V? = = F 
a 


where M is the particle feed rate; Af = 


L 
v’ and U is the 
avcrage velocity of flow. 


where z is the percentage of i greater than and 
Vue Vu 

the probable value of z is given in Fig. 4, obtained from 

statistical velocity distribution developed from kinetic 

theory of gases [14). 

Ha.istTrom [15] measured experimentally the maximum 
carrying capacity without build-up of a turbulent stream 
in a 3in. diameter duct at various average stream velo- 
cities, using glass beads of sizes from 105-123 microns in 
diameter. The intensity of the stream can be calculated 
from the above equation with the aid of Fig. 4: 
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Fic. 4. Probable velocity distribution. 





Measured 
Quantities Calculated Quantities 

Mean dis- 
Air flou 


rate, 


tance 
U M 
[ps lbs /min 


Particles’ between 


3 


Vu2/U 
lbs / min in particles, 
dia- 
meters 

— oa 
262 | O-T27 
28-8 1-26 
34°8 8-12 





APPENDIX II 
Integration Formula Developed for the Approximate Solution 
oO @ 
¢ a2n2 ‘ 
- dn 


Int ¢ 
: J aind;) b4*n? j b*n* b%n8 b8 8. \dn 
+ bn 


2a 2a? 
. 9 
= G (b? a”) ‘Van 2a 
b?n? 


~ cos 2 pn dn b*n* 


)cos 2 pn dn 
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The general form of the above integral is 


@ 
con f . » 
—~____ cos 2 pn dn = Te r'le 
1 + b'n? 2a 
0 


10” 180° 
1- + - 
| 2a* 27 a* 
“(5 820 8523 0° 8-5-7 2-8-4 b* 
a\a* 2ia® 2120 2 128° 
p*/o* sad 57340° 5-793-450'° 
An 210° 2120 2 1-280" 


1-3-5 b® 
2° a® 


Na 2ia 2 120° 2 i2sa? 
p?/b® 956% 911565" 9-11-18 56-7 5"* 
@\ 21a ie? 2 
=(5-2 6b 1 13 6-7 b'¢ ) : 
10 


+- . 
a’®\a 21a 2° 12a'* 


Vere n (5 5) 


p?(o® 740 70450" 7-9-1145-6b"? 
| 

1-2-8 al4 
“_ 


J 


2a a a 


= V"—pi/e | (5) + mi (5 


2a a 


NOTATIONS 


(with dimensions of each quantity in mass, M ; length, L ; 
and time, 7) 


A amplitude of the nth harmonic of velocity fluc- 
tuation, LT" 
Cp = drag coefficient, dimensioniess 
d, = diameter of particle, L 
E = eddy diffusivity, L*?T~' 
E, = diffusivity of the particle, Lr 
F = 18 »(d,* py), T" 
J(n) = Eulerian statistical frequency function, T 
J, (mn) = Lagrangian statistical frequency function, T 
Fy = root-mean-squared Froude number +/wu?/d,g, 
dimensionless 
G(K) = function of K defined by eq. (29), dimensionless 
g = Gravitational acceleration, LT* 


1(K) = function of K defined by eq. (81a), dimensionless 
K = defined by eq. (80), dimensionless 
k = an empirical constant 
1, = Lagrangian scale of turbulence of the stream, L 
,’ = Lagrangian scale of turbulence of the particle, L 
l, = Eulerian scale of turbulence, L 
m <= mass of a particle, M 
= frequency, Tt 
R, = Eulerian correlation of the stream motion, dimen- 
sionless 
Lagrangian correlation of the stream motion, 
dimensionless 
Lagrangian correlation of the stream motion, 
based on », dimensionless 
‘ = Lagrangian correlation of the particle motion, 
based on » or »’, dimensionless 
Reynolds number of particle, dimensionless 
dy up 
root-mean-squared Reynolds number, -” 
dimensionless 4 
long time interval, T 
time, T 
= axial component of turbulent velocity component 
of the stream, LT=' 
vertical component of turbulent velocity of the 
stream, LT" 
= Axial component of turbulent velocity of the 
particle, LT! 
vertical component of turbulent velocity of the 
particle, LT! 
relative velocity between the particle and the 
stream, LT" 
= intensity of turbulence of the stream, L?7-* 
= axial intensity of turbulence of the particle, L? 7? 


urs 
= intensity of the particle relative to the stream, 
Urs 
p = Viscous dissipation of the particle, per unit time, 
ML*T* 
= measured in the axial direction of flow, L 
= average displacement, L? 
= angle of phase lag of velocity of particle from the 
stream, dimensionless 
» = defined by dy = 4/u? dt, L 
9’ = defined by dy’ = 4/u,? dt, L 
A = defined in reference [7] 
A,, = defined by eq. (7), @ measure of size of eddy of the 
stream, L 
A’, = @ quantity analogous to A,, but for the particle, L 
pp = viscosity of the stream, ML~'T-! 
p = density of the stream, ML~* 
Pp = density of the particle, ML* 
¢ = time measure, T 


(Note—Notations in the Appendices are explained in the 


Appendices.) 





(1) 


(2) 


[3] 


Statistical properties of momentum transfer in two-phase flow 


REFERENCES 


Nicnoison, H. M. and Fiexp, J. P.; Some experimental techniques for the investigation of the mechanism of 
flame stabilization in the wakes of bluff bodies. Third Symposium on Combustion, Flame, and Explosive Pheno- 
mena 1949 45. 


LATINEN, G. A.; Mechanisms of fluid-phase mixing in fixed and fluidized beds of uniformly sized spherical par- 
ticles. Ph.D. thesis in Chemical Engineering, Princeton University, 1951. 

Lonowe.i, J. P.; High Speed Aerodynamics and Jet Propulsion, Vol. 2, Combustion Processes, Princeton Uni- 
versity Press,1956. 

Vanont, V. A. and associates; Dynamics of particulate matter in fluid suspension. Report No. N-71.1, Hydro- 
dynamic Laboratory, California Institute of Technology, November, 1950. 

Pinkus, O.; Pressure drops in pneumatic conveyance of solids, J. Appl. Mech. 1952, December, 425. 

Lin, C. C.; On the motion of a pendulum in a turbulent fluid, Quart. Appl. Math. 1943 43. 

Taywor, G. I.; The spectrum of turbulence, Proc. Roy. Soc. (London), 1937 164 426. 


Taywor, G. 1.; Statistical theory of turbulence, Proc. Roy. Soc. (London) 1935 151 421. 


Tayvor, G. I.; Diffusion by continuous movements, Proceedings of the London Mathematical Society 1922 20 196. 


Bassett, A. B.; A Treatise on Ilydrodynamics, 1888 George Bell & Sons. 

Gripert, M., Davis, L. and Atrman, D.; Velocity lag of particles in linearly accelerated combustion gases. 
Progress Report No. 20-194, ORDCIT project, Jet Propulsion Laboratory, California Institute of Technology, 
May 28, 1953. 

Soo, S. L.; Wet-compression in an axial-flow compressor, Trans. A.S.M.E. 1952 74 879. 

Drypben, H. L.; A review of the statistical theory of turbulence, Quart. Appl. Math. 1943 1 7-42. 
Kennarp, E. H.; Kinetic Theory of Gases, 1938, McGraw-Hill Book Co., Inc., New York, 51. 


Hatistrom, E. A. N.; Design of experimental apparatus for the study of two-phase flow in circular straight 
pipe. M.S. thesis in Mechanical Engineering, 1953, Princeton University, Princeton, N.J. 





Chemical Engineering Science, 1956, Vol. 5, pp. 68 to 80. Pergamon Press Lid. 


Heat transfer to boiling binary liquid mixtures 


W. R. van Wuix, A. S. Vos and S. J. D. van STRALEN 


Laboratory of Physics and Meteorology, Agricultural University, Wageningen, The Netherlands 


Abstract—Heat transfer in boiling liquid mixtures of water and methylethylketone has been 
measured throughout the ranges of convection and nucleate boiling at atmospheric pressure [43] 
With some mixtures a consi. erubly higher maximum heat flux in nucleate boiling has been 
found than in the case of water. This higher heat flux is obtained at the same temperature of 
the heating surface as for water, or alternatively, the same heat flux is obtained at a lower 
surface temperature. 


For mixtures of water with acetone, methylethylketone, the lower aliphatic alcohols and 
ethyleneglycol respectively, the maximum heat flux in nucleate boiling has been studied as a 
function of concentration. In addition, mixtures consisting of two organic Iqiuids, i.e. dioxane- 
methanol and 2-chloroethanol—di-iso-propylether, were investigated. In all mixtures a maximum 
value of the maximum heat flux for nucleate boiling occurs at a certain concentration. The higher 
the number of carbon atoms in the alcohol and ketone series, the more this maximum is shifted 
towards lower concentrations of the organic component. A maximum of 3} times the value 
for water has been observed in the s) stem water—1-pentano]. Partial miscibility has no influence 
on the continuity of the maximum heat transfer curve. 


The occurrence of the maximum is qualitatively explained by the higher dew point of the 
vapour formed by flash vaporization, as compared with the boiling point of the original liqiud 
mixture. At a certain concentration, the dew point of the vapour bubbles muy become equal 
to the temperature of the surrounding superheated liquid, and in that case the vapour bubbles 
grow only by diffusional mass transfer The average size of the bubbles leaving the heating 
surface is then smaller, and the maximum heat transfer is higher. For mixtures of water and 
alcohols, it has been calculated that 1 mole % 
surface. For water—acetone and water—methylethylketone, only 0-3 and 0-1 mole °® 


of the liquid is immediately vaporized at the 
direct 


v 


vaporization respectively is found 


Résumé— Les auteurs ont étudié et mesuré le transfert de chaleur dans des mélanges a !'ébullition 
d'eau et de methyléthylcétone (43) d'un bout a l'autre de l'intervalle convection — nucleate 
boiling. 


Pour certains mélanges ils ont trouvé un maximum du flux de chaleur dans le ‘‘nucleate boiling” 
considérablement plus é¢levé que dans le cas de l'eau. Ce flux de chaleur plus élevé est obtenu 
pour la méme température de surface chauffante que dans le cas de l'eau ou réciproquement le 
méme flux de chaleur s‘obtient avec une température de surface plus faible. 


Dans le cas de mélanges d'eau et d'ecétone, de methyléthyleétone, d'éthyléne glycol et d'alcools 
aliphatiques inférieurs respectivement, les auteurs ont étudié le maximum du flux de chaleur 
dans le * nucleate boiling *' en tant que fonction de la concentration. En outre, ils ont experi- 
menté des mélanges de deux liquids organiques, par exemp1> dioxane — methanol et 2-chloroethanol- 
di-iso-propylethyléther. Pour certaines concentrations, le flux de chaleur maximum pour 
le “ nucleate boiling " présente un maximum. Plus le nombre d’atomes de C est grand dans 
les alcools et cétones, plus le maximum est déplacé vers les concentrations faibles en composant 
organique. La valeur de ce maximum atteint trois fois et demi la valeur observée avec l'eau, 
dans le cas du systéme eau-pentanol. La miscibilité partielle est sans influence sur la continuité 
de la courbe du transfert de chaleur maximum. 


La production de ce maximum s'explique qualitativement par le fait que le point de rosée 
de la vapeur formée par vaporisation instantanée est plus élevé que le point d'ébullition du 
mélange liquide primitif. 


A une certaine concentration le point de rosée des bulles de vapeur peut devenir égal a la 
température du liquide surchauffé environnant et dans ce cas, les bulles de vapeur ne croissent 
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Heat transfer to boiling binary liquid mixtures 


que par transfert de masse par diffusion. 


chauffante est alors plus petite et le transfert maximum de chaleur est plus élevé 
mélanges d'eau et d'alcools, les auteurs ont calculé qu'une mole ° 


La dimension moyenne des bulles quittant la surface 


Pour des 


de liquide est immédiatement 


° 


vaporisée a la surface, pour les mélanges eau-acétone et eau-méthyléthyleétone ils n'ont trouvé 


que 0,3 et 0,1°% moles directement vaporisées 


GENERAL SURVEY 


Heat transfer to boiling liquids has been the 
subject of many investigations. It is discussed 
in current textbooks, e.g. of McApams [83] and 


of Jako [20]. 


—_ 


Micro 
ton COrwection 


Filmbor ng- 
Radiation k 
/ 
/ 


Partial 
f imboiing 


Nucleate 
boiling 


comemmns 


. 


o— 


| fy 
Fic. 1. Schematic representation of heat flux q/A as a 
Junction of temperature difference 0 between heating surface 
and bulk of liquid. 


The curve showing heat flux ¢/A as a function 
of temperature difference @ between heating 
surface and bulk of liquid is schematically given 
in Fig. 1. 
first of which is OA in which heat is transferred 


Five regions can be distinguished, the 


by convection (and also by conduction), and no 
vapour bubbles are observed. This also holds 
true for the region AB. The convection in the 
latter region is, however, much more intense than 
in the region OA. We are, therefore, inclined to 
speak of 
as distinct from 
point B formation of bubbles starts on special 


Tr 


micro "’-convection in the region OA 
“macro "’-convection in AB. At 


spots of the heating surface, and nucleate boiling 
takes place in the region BC. The number of 
nuclei at which vapour bubbles are generated, 
and thus the fraction of the surface covered with 
vapour bubbles, increases in the region B to C, 
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until in the region CD a gradually increasing area 
of the heating surface is covered with a coherent 
At the same time the heat flux 


drops as a result of the presence of the insulating 


layer of vapour. 


vapour layer, although the temperature difference 
@ between heating wire and bulk liquid increases. 
Thereafter the heat flux rises again, owing to the 
until the 


contribution of radiation, 


E. 


increasing 


wire burns out at 


Nucleate boiling heat transfer has been investigated 
extensively and mostly for water. The maximum at C 
1) has been mentioned by several authors, mainly 


11, 28, 24, 34, 37, 41] 


(Fig 


for water (4, - 


Comparatively little work has becn done on mixtures. 
Bonitta and Perry [4] have investigated some binary 


mixtures of water with ethanol, acetone and butanol 


respectively, as well as binary mixtures of ethanol—butanol 


and ethanol—acetone, CrcHeLu and BoniiLa [7] have 


investigated mixtures of water—ethanol and propane—n- 
5] water—styrene and 


[13] 


ethanol, benzene 


BoniL_a and EISENBERG [ 
butadiene ; Bonet1 

furfural ; Kirschbaum water- 
toluene [28] and water—glycerol [29] 


pentane ; 


water and GESTER hydro- 


carbons 


In a narrow concentration range in which water was 
present in excess, experiments were made by CrYDER 
and GILLILAND [9] on glycerol, saccharose, NaCl, NagSO, ; 
on glycerol, NagSO, ; by 
by Ruopes and BripGes 


) 


by Cryper and FINALBOoRGO [8] 
LARSON [30] on NaOH, CaCO, ; 
[41] on oleic acid with NaCl and NagCO,; by InsINGER 
and Buss [18], Jakos and Linke [24], KirscHBaumM 
[29] and Morcan ef al. [36] on wetting agents; by 
McApams et al. [35] on air. 


In some of their mixtures BonrLua and Perry have 
found a maximum heat flux in nucleate boiling (at point C 
in Fig. 1) exceeding somewhat that of either one of the 
pure components [4]. No systematic investigation about 
the influence of concentration was, however, made, and 
the increase of maximum heat flux mentioned by them 


was very moderate. 


Kirscupaum [29] has found that in 20% by weight 
solution of glycerol in water the overall coefficient of heat 
transfer was twice as high as with pure water, or as in 
a 50% by weight solution of glycerol, at the same tem- 
perature difference 9 = 20°C. He attributes this behaviour 
to foaming. No sufficient data are, however, given to 
conclude that the maximum heat flux was reached in 


this case. 
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Table 1. Survey of data at atmospheric pressure concernng pow 
Com ponent Bowling point Azeotrope Solén A 
%~B 
A B T,4 Tr TAs —_—— _— T A-T Az\ Reference 4/100g Tem 
| wt Mole Wee 
| | | -—<s 
water | ethanol 100 78-8 78-15 05-57 80-41 | 21-85 17, 27, 81 ~ ~ 
78-8 78-18 95 59 89-45 21-87 81 
o 1-propanol 100 97-2 87-72 71-60 43-17 12-28 17 ~~ ~ 
| 97-19 87 72 71-69 43-17 i223 6] 
| ora 87-75 71-95 43-40 12-25 81 
———EEE - aa = 
2 | butanol 100 | 82-55 79-01 88-24 64.60 20-09 17, 81 ~ - 
79-96 90-8 70-6 81 
™ 1-butanol 100 117-75 92-4 62 28-4 7-6 17 7-9 bu 
| 117 92-4 62 23-4 7-6 7,81 7-8! a 
| | 1160 02-25 63 29.8 7-75 81 7 33 cal 
| 117-70 | 92-50 57-5 24-76 7 50 26, 27 “ a 
| 7” v2 
S76 | 68 40, 42 ) im 
62-67 2 1) a 
— = —— = — 
1 -pentanol 100 187-8 95.05 ae 148 4-05 17 27 26 oa 
| 187-8 05 05 466 15-1 4-05 81 2-26 2 
187-7 05-4 “6 148 46 81 
| 188-0 96.0 46.0 148 40 81 
- 1-hexanol 100 157-85 7-8 25 5 54 2-2 17 0-628 5 2 
j | 0-58 54 at 
—— - — - oe 
” 1-heptanol 100 176-15 08.7 17 808 1-8 17 0-280 -oee4 | 10 
| 0-180 @-1 276 2 
| 0-08 0-o1l# 1! 
” l-octanol 100 105-15 90-4 10 1-5 0-6 17, 31 00-0538 0-00745 2 
|} 195-15 90-5 11 1-7 0-5 81 
$$ _________ fee So —— . — - — — : 
a | acetone 100 ae non-azeotrope | 17, 81 ~ = 
56 50 non-azeotrope 27 
= ° _ Lo . ae aa — - } 
' ’ 
- methylethylketone 100 | 76 73-45 89 660 26 55 17 19 15 of 9 
70-6 73-45 88-6 66.0 26-55 81 15-76 S 
72-56 | 78-57 88 62 6607 | 26-48 81 18-0 7 
760 | 736 89-1 e7-1 86| | (264 a7 22 of? 2 
| 70-6 73-8 7 aa 85-8 26-0}! 
726 86| «(7386 | 24 | 82 26 sf 2 
| | 89 | oo | :. - as-ai6 ‘ 
j | wo 16 ¥ 
| | a7-s : 
j ' | 
= = ———— = — cers — 
dioxane methanol 101-4 647 non-azeotrope 17 — " 
101-05 64.60 non-azeotrope 27 
_— aa -— | + Samaria ES — 
2-chloroethanol | di-iso-propylether 128 68-5 non-azeotrope | 17 ~~ ~ 
| 128-00 68-40 non-azeotrope 27 
j i — 
ethyleneglycol | water (747 mm)| 196-70 99.60 non-azeotrope 27 ~ o~f 
| 197-85 100 non-azeotrope 17 
197-4 100 non-azeotrope | 17,81 
a — - —— _——_— _ —— ’ — — 
1-butanol water | 
1 pentanol water | 
i be 
In systems for which no complete equilibrium data are known, the equilibrium constant has been estimajm the 
its ratio to the value for component A are given. The experimental dapfis ma 
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g point, azeotrope, solubility and maximum heat flua. 














4 Heat flux (q/ A) maz Data at maz. (q/A) max Data at theor. maz for G selected 
. . ee . . = — meen 
Maximum % B % B Selected 
Temp Ref K Maz A Ratio T K — . T K SS G 
wt Mole Wi Mole 
—— | 
10-2 80 18 1:7 86 48 | 20 91-7 7-2 10 + 0-01 
= — — Ee =_ 2 — 2 EE 
11-6 34 17-6 1:9 89-2 51 18 7 93 10-7 6 2 0-01 
| 
= 7 = cibeusenstiaill — —— EE 
43 20 22 24 0-6 
ve 19 24-7 45 17-6 2:5 07-25 21-5 1-6 0-4 05-5 19-5 3 0-75 0-01 
a 15, 19 (> 15) 
row 16,40 
25 42 
4 15 
La! 19 
a 15,19 
22 16 (> 82) 55 16 34 1-8 0-35 
25 ra 
25 «0 ( 53) | 
20 16 
‘ = a : — 
100 16,40 (> 120) 
25 40 | 
18 16 
_ «0 (> 205) 23 13-5 1-7 (99-7) 0-35 0-05 
= — 3 - EEE . EE 
~ 28 28 12-8 1:8 78 145 12 4 84.8 24 7 2-3 0-003 
| 
. | ———— Se EE —E 
90 16,40 +5 41 17-7 23 89 87-3 4 10 86 bad 5 1-8 0-001 
80 | 19 (> 15) | 
736 | 382 | 
20 19 
10 16,40 
20 } ae | 
80 | 19 | | 
2 | | 
20 89 
— -——_ — —_———__ >—_——_- —_— - — — ~ — — + — _—— 
14 15-6 i124 13 75 9-7 2 5 75 9-7 2 5 
89 22-6 11-7 19 90-7 11-6 s 65 110-5 81 2 1-5 
10 21-7 17 13 174 y 2 7 171 8-3 2-5 ~ 
_— = — — -—- — _ — —- + -- —_}— - — = — }—— -— -— —— —_—— = 
16-5 18 11 12 100 1-7 6 21 111-4 7 1 4 
—_ ——_—— S — — —___—— _—————<——| —___ wy $$} $$$ }_— —EEEEE 
11 8-5 13 | @ 16-5 




















the relation K = y4,/z,,.. The maximum of maximum heat flux in nucleate boiling for mixtures, and 
maximum are compared with their values predicted by theory. 
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W. R. van Wok 

In the present article the authors have system 
atically investigated the influence of concentration 
upon the maximum heat flux in nucleate boiling 
Vos and vAN STRALEN [43] had previously found 
that for some mixtures of water and methy! 
ethylketone a maximum heat flux occurs equal 
to twice that of water 

The present experiments were carried out at 
atmospheric pressure with the same boiling vessel 
as described previously [43], but the horizontal 
diameter 0-02 cm 


platinum wires (purity 99-99% 
and length about 6 cm) were heated with alter- 
The difference 


between wire and bulk liquid was not measured 


nating current. temperature 


Only the voltage and the current through the 


heating wire were measured in order to determine 


the rate of heat transfer q. 


Mixtures of water with acetone and methyl- 


ethylketone, of water with ethanol, 1-propanol, 


t-butanol, 1-butanol, 1-pentanol and l-octanol 


respectively, and of water with ethyleneglycol 


have been investigated. Experiments were also 


made with binary mixtures consisting of two 


» 


organic liquids, i.e. dioxane—methar and 2 


chloroethanol—di-iso-propylether 

The 
and were 
chemically the the 


1-butanol had a boiling range of 115 to 117°C, 


used in all 
the 


latter furfural! 


Distilled water was 
methylethylketone 


cases 
1-pentanol 


pure, free, 


and the other liquids were of commercia! grade 


The data of mutual solubility at room temperature 

were taken from Internationa! Critical Tables [19)}, from 
The solubilities at 
the azeotropic boiling point are only known for mixtures 
of water—methylethylketone [16, 19, 82, 40) and water—1 
butanol [15, 19) 
mixtures were taken [ 
further: ethanol (25), 1-butanol [42] 
{8, 32, 88, 29). Data pertaining to t-butanol 
and l-octanol have not found Azeotropi 
have also been taken from references [17, 31 
of these data is shown in Table ! 


HopoMan [16] and from Peary [40) 


The boiling and equilibrium data of the 
from literature [26, 27, 40) and 
methylethyiketone 
1-pentanol 


data 
A summary 


been 


RESULTS 


Miztures of water 
acetone 


methylethylketone and water 


been 
another paper [48]. The curves showing the heat 
flux g/A as a function of the temperature differ 


Some measurements have described in 


4. S. Vos and S. J. D. 


VAN STRALEN 


ence @ between heating surface and bulk liquid 


for boiling mixtures of water and methylethy]l- 


ketone are represented in Fig. 2. The mixtures 


with 4-2 and 88-5% by weight methylethylketone 


‘ 0 x00 
homogeneous, those with 20% and 52%, are 


are v 


in the heterogeneous region. 


—— VYAlcal sec'cm 


s 


—= 6 (XL 


Fic. 2. Water—methylcthylketone. Heat flux q/A for con- 

veetion and nucleate boiling as a function of temperature 

difference @ between heating surface and liquid. The 

figures near the curves denote % by weight methylethy! 
ketone. 


below @ 


at a higher temperature differ 


Micro-convection occurs Lo" 
approximately ; 
ence a stronger convection sets in until the first 
visible bubbles occur at @ = 25 
The 
tinued until filmboiling started. 
the wire burned out. The size of the bubbles 
leaving the heating surface is definitely smaller 
for 42% methylethylketone than with the other 


muxtures or with the pure components. 


‘ apour 


approximately measurements were c 


In some cases 


—= 10% YA (piu prin? 
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The maximum heat flux (at point C in Fig. 1), 
which forms the separation between the regions 
of nucleate boiling and filmboiling, differs con- 
different 
highest for the mixture containing 4-2% methy]- 


siderably for concentrations. It is 


ethylketone and it is then 150% higher than for 
pure water. 


7 


—=— GA max ical sec em? 
; 8 
| 
| 
> 


\& 
r) 
> 
E 
uu 
“ 
o 
AZ 
3 
< 


“130 
_——- *. @-7 s 


Fic. 3. Water 
and maximum coefficient of heat transfer h = ¢ 


methylethylketone. Maximum heat flux ¢/ A 
A @ for 
nucleate boiling as a function of temperature { = 6 + T 
of heating surface. The figures near the circles on the 
The 


dotted lines serve for a comparison of the mixtures with 


curves denote % by weight methylethylketone 


the pure components 


The maximum heat flux g/A and the maximum 
coefficient of heat transfer h = ¢/A @ are repre- 
sented in Fig. 8 as a function of the temperature 
T, where T denotes 
The 


maximum heat flux which is 24 times the value 


of the heating wire t = @ 
the boiling temperature of the liquid. 


for water occurs approximately at the same 
temperature of the wire as with pure water. In 
another mixture containing approximately 50% 
methylethylketone the maximum heat flux is 
equal to that of water but the temperature of the 


since tre 


73-U. 


wire is approximately 20°C lower 


boiling temperature 7 is approximately 


This mixture consists of two yuid phase Ss 


—e Mole M Kefone 


6U 





—e Wi %/> ME Ketone 


Fic. 4. Water 


as a function of composition 


methylethylketone. Maximum heat flux q/ A 
Measurements carried out 
with the same wire are represented by the same figures 
The dotted 


boundaries of the region 


Az azeotrope vertical lines indicate the 


of demixing at azeotropic boiling 


point. The other curves represent AT as a function of 
composition for a constant vaporized molar fraction G 
\T is the of the 


vapour bubbles and boiling temperature of the original 


difference between dew temperature 


liquid. The numbers near these curves are the values of G 


Maximum heat flux as 


composition 1s 


a function of liquid 
Fig 4 The dotted 
boundaries of the 


It is interesting to observe 


shown it 


vertical lines denote the 
region of demixing. 
that no discontinuity is found in maximum heat 
data 


Ww here 


transfer when demixing occurs Some 


on GSC shown in Table 1, 
with 


concerning (¢ A) 
7 


a comparison other binary mixtures is 
also made. 

Corresponding measurements for the system 
A less 


pronounced maximum appears at about 12% by 


water—acetone are represented in Fig. 5. 


weight acetone. A second maximum is observed 
at a higher concentration. This maximum cannot 
be explained by the theory which will be given 


below. 
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Fic. 6. Water—ethanol. Maximum heat flux ¢/A as a func- 

tion of composition. Measurements carried out with the 

same wire are represented by the same figures. Az azeo- 

trope. The other curves represent AT as a function of 

composition for a constant vaporized molar fraction G. Fic. 8. Water—t-butanol. Maximum heat flux ¢/A as a 
AT is the difference between the dew temperature of the function of composition. Measurements carried out with 
vapour bubbles and boiling temperature of the original the same wire are represented by the same figures. Az = 
liquid. The numbers near these curves are the values of G. azeotrope. 
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Heat transfer to boiling binary liquid mixtures 


Miztures of water with some lower aliphatic alcohols —= Mole % 1-Pentanol 


Maximum heat flux for nucleate boiling was also 


determined as a function of concentration 
ethanol (Fig. 6), water—1- 
t-butanol (Fig. 8), 


1-pentanol 


mixtures of water 
propanol (Fig. 7), water 
water—l-butanol (Fig. 9), 
(Fig. 10) and water—l-octanol (Fig. 11). 
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—e Wi Y% /- Butanol 


Water-1-butanol. Maximum heat flux ¢g/A as a 
function of composition. Measurements carried out with 
the same wire are represented by the same figures. Az 
azeotrope. The dotted vertical lines indicate the boundaries 
of the region of demixing at azeotropic boiling point. The 
other curves represent AT as a function of composition 
for a constant vaporized molar fraction G. AT is the 
difference between dew temperature of the vapour bubbles 
and boiling temperature of the original liquid. The num- 
bers near these curves are the values of G. 


Fic. 9. 


Increasing the number of carbon atoms in the 
alcohol series shifts the maximum towards lower 
alcohol concentrations. This maximum becomes 
higher and steeper from ethanol up to 1-pentanol, 
but it then decreases and for 1-octanol it is back 
again at the value of ethanol (compare Table 1). 

With 1-pentanol the very high maximum value of 
55 cal sec“ cm™* at 1.8% by weight pentanol is 
obtained at a boiling point of about 96°C. This is 
about 8} times as high as the maximum for water. 
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10. Water—1-pentanol. Maximum heat flux g/A as a 
Measurements carried out with 
the same wire are represented by the same figures. Az = 
azeotrope. 
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Fic. 11. Water 
function of composition. 
the same wire are represented by the same figures. 


1-octanol, Maximum heat flux ¢/A as a 
Measurements carried out with 


With 1-butanol and second 
maximum at high alcohol concentrations appears. 
Again the curves do not show a discontinuity 


or a break if a second phase and emulsification 


l-pentanol a 
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appear. This might be explained by assuming a 
preferential wetting of the heating surface by 
one of the phases. 

It may be of importance to stress the fact that 
in all mixtures described in this paper no stable 
however, at certain 
of the 


definitely smaller, which causes more 


foaming was observed ; 


ranges of mixtures the size bubbles is 


tendency 
to foaming. 
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Fic. 12. Correlating maximum nucleate boiling heat fluc. 

(a) Correlation by Cicne.ti1 and Boni.1a ; 

ethanol (the figures denote % by ethanol) ; 
(c) water—1-pentanol (1:8% by weight 1-pentanol). 


(6) water 


weight 


With the data of critical pressures p. for water—ethanol, 
determined by Griswoitp [14] we plotted in Fig. 12 
(G/A)max/Pe #8 & function of reduced pressure p, = p/p, 
in the same way as has been done by Cicwe.tu and 
Boniiia [7]. The dotted line (a) represents their corre- 
lation for several pure components and some mixtures 
at various pressures [4, 7]. This curve (a) shows a maximum 
at p, ~ 03, which is not shown in Fig. 12. 
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The curve (b) — where the figures denote % by weight 
our data 


deviates from the curve (a) but shows essentially the same 


ethanol - represents for water—ethanol, and 
shape as the curve in Fig. 6; the maximum, however, is 


shifted now to a higher concentration 


For a 18% by weight mixture of 1-pentanol in water a 
(¢/A) may Of 55 cal sec”! cm~* was observed. Although no 
critical pressures for this system are found in literature, 
an estimation of the situation of that mixture in Fig. 12 
can be made, as it is reasonable to assume that the critical 
pressure for this mixture does not exceed that of water. 
This mixture namely shall be situated then on that part 
(c) of the straight line through the origin with a slope of 
. P/Pe 
This means that for 


55 cal sec"! cm~* atm, that is characterized by 
> 0-0045 (the value for water). 
this mixture (¢/A),.6./P¢ 


as in the correlation. 


might be at least twice as high 


In the mixtures discussed the most pronounced 
maximum was always found to occur for a con- 
centration in which water was in excess. A 
qualitative explanation of this phenomenon is 
the 


volatility of the constituents in the mixture. On 


given below on basis of the difference in 
the strength of this theory a second maximum is 
expected to occur for a low water concentration, 
if water is considerably more volatile than the 
organic compound. The position of the second 
maximum for a low water concentration is in 
agreement with theory in case of 1-butanol and 
1-pentanol. 

If the mixture forms an azeotrope at a certain 
concentration, water is the more volatile com- 


If no 


azeotropic concentration exists, water is either 


ponent on one side of the azeotrope. 


the less volatile or the most volatile component 
throughout the entire range of concentrations. 


The latter situation is realized in mixtures of 
water and ethyleneglycol. 


Mixtures of water—ethyleneglycol 

In accordance with the theory only the second 
maximum heat flux when ethyleneglycol is in 
excess, is observed. This is shown in Fig. 18. 


Miztures of two organic components 


Two binary systems (dioxane—methanol and 
2-chloroethanol—di-iso-propylether) each con- 
sisting of organic compounds only, have been 
investigated to ascertain that the observed 
increased heat flux is not due to some property 
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as a function of composition. Measurements carried out 
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difference between dew temperature of the vapour bubbles 
and boiling temperature of .the original liquid. The 
numbers near these curves are the values of G. 
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—- Wt Yo Methanol 

Fic. 14. Diorane—methanol. Maximum heat flux ¢/A as a 
function of composition. Measurements carried out with 
the same wire are represented by the same figures. The 
other curves represent AT as a function of composition 
for a constant vaporized molar fraction G. AT is the 
difference between dew temperature of the vapour 
bubbles and boiling temperature of the original liquid. 
The numbers near these curves are the values of G. 


characteristic for systems that contain water as 
a constituent. Here again qualitative agreement 
In Fig. 14 the 


curve for dioxane—methanol shows a little pro- 


with theory is actually obtained. 


nounced maximum at an excess concentration of 
An 


was 


additional maximum, however, 
the 


cannot be 


dioxane. 


appears—as also case in water 


acetone mixtures—which explained 
by difference in volatility of the constituents in 


for 
Fig. 15 


the mixture. The curve 2-chloroethanol 


di-iso-propylether in shows a more 


pronounced maximum. 


—e Mole % di-1-Propylether 
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Fic. 15. 2-Chloroethanol—di-iso-propylether. Maximum heat 
flux ¢g/A as a function of composition. Measurements 
carried. out with the same wire are represented by the 
same figures. The other curves represent AT as a function 
of composition for a constant vaporized molar fraction G. 
AT is the difference between dew temperature of the 
vapour bubbles and boiling temperature of the original 
liquid. The numbers near these curves are the values of G. 


QUALITATIVE EXPLANATION OF 
MAXIMUM HEAT TRANSFER 


It has been shown by other investigators that 
most of the heat which is ultimately consumed in 
vaporization of water is first transferred to the 
bulk liquid in the neighbourhood of the heating 
surface, and only a small] fraction is used for a 
direct formation of bubbles in contact with the 
heating surface [21]. The bulk liquid is slightly 


superheated and the small vapour bubbles 


77 





W. R. van Wu, A. S. Vos and S. J. D. 


released from the heating surface grow rapidly 
during their passage through the liquid. Film- 
boiling starts if, owing to a larger temperature 
difference @ between heating surface and liquid, 
a vapour film is formed at the heating surface 
if the vapour 
bubbles merge in the neighbourhood of the 


(Leidenfrost phenomenon), or 
surface. Thus the smaller the average size of the 
bubbles leaving the heating surface, or the better 
the wetting, the more the onset of filmboiling is 
shifted towards higher heat flux. 

If only one component is present, the vapour 
bubbles, which are assumed to be at their dew 
point, must always grow during passage through 
a superheated liquid; but this is no longer so in 
case of a mixture. 


Let y and z be the molar fraction of the most 
volatile component in the binary mixture in the 
vapour phase and in the remaining liquid after 
Consider a system in 
the 
organic component is more volatile than water. 


vaporization respectively. 


which in a region of low concentrations, 


For equilibrium flash vaporization one has : 
y=Kze (1) 


in which K is the equilibrium constant of the 
most volatile component. Let further z, be the 
molar fraction in the liquid phase before vapori- 
zation has taken place. Then, if a molar fraction 


G is vaporized, one has the relation : 
(l—G)z+Gy=z2, (2) 
and from equations (1) and (2) : 


r=2,/1+(K —1)Gand y=Kz,/1+(K—1)G (8) 


The dew temperature of the vapour bubbles is 
equal to the boiling temperature of the remaining 
liquid and since z is smaller than z,(as K > 1) 
for the volatile component, the dew temperature 
of the bubbles 7 + AT exceeds the boiling 
temperature 7 of the original liquid. The dew 
temperature depends on the value of G. It can 
be calculated from equation (8). It is nearly 
equal to the boiling temperature of the original 
liquid if G is small and if G = 1 it is the tem- 
perature at which a vapour containing a molar 
fraction z, condenses. 
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Thus if the bubbles rise through a mass of 
superheated liquid of the original composition, it 
may happen that a certain value of G the dew 
temperature equals the temperature of the liquid. 

One has then a situation in which the vapour 
in a bubble is in equilibrium with a liquid bound- 


ary layer which possesses a composition different 
from the surrounding liquid, but which has the 
same temperature as the latter. No heat flows 
from the superheated liquid to the boundary 
layer; and it is only owing to diffusional inter- 
change between the boundary layer and the 
surrounding liquid that composition and dew 
temperature of the vapour will change and 
further vaporization can occur. A smaller size 
of bubbles resulting in a larger maximum heat 
transfer may then be expected. 

Along this line of thought the occurrence of 
most of the maxima in the binary mixtures can 
be explained. In all cases investigated the size 
of the bubbles was actually a minimum at the 
As an example, the size 
1-butanol 


maximum heat flux. 
of the bubbles is shown for water 
in the photographs of Fig. 16. 

The dew temperature of the vapour minus the 
boiling temperature of the original liquid AT is 
plotted as a function of z, in Figs. 4, 5, 6, 7, 9, 18, 
14 and 15. The value of G is constant on each 
curve. This dew temperature and the boiling 
temperature of the original liquid were read from 
graphs of the boiling temperature of the binary 
mixture as a function of z. The corresponding 
values of K were calculated from a curve giving 
y as a function of z. 

A maximum value of AT occurs in all curves. 
In the case of the alcohols, these maxima occur 
at the same value of z, as the maximum heat flux 
The difference AT 
is then some tenths of a degree centigrade. This 


if G is approximately 0.01. 


is the same order of magnitude as found by other 
investigators for superheated pure water [22]. 
Since mixtures in the neighbourhood of the maxi- 
mum heat flux consist mainly of water this is a 
satisfactory result. In mixtures of water—methy]- 
ethylketone a maximum value of AT as a function 
of concentration is found at a value of z = 0-015 
(approximately 5% by weight). The position of 
this maximum is hardly dependent upon the 


























ic. 16. Water 1-butanol. Photographs showing the size 


of vapour bubbles at constant heat flux (7-5 cal sec lem) 


for: 1, pure water, 2. water with 2% by weight 1-butanol, 


ind 3. pure 1-butanol. The length of the heating wire is 


Gem and the diameter 0-02 cm. 











Heat transfer to boiling 


alue of the vaporized fraction G. It is, therefore, 
the 
by comparing the position of the 
calculated value of AT with the 
If, how 


some tenths 


mpossible to select a value of G as most 
probable one, 
maximum of the 

tual position of maximum heat flux. 


\T ma 


; 


ximum is assumed to be 


gree centigrade, a value of G 0-001 Is 


ilculated. This is ten times smaller than in the 
se of the 


0.003 is determined and for the other systems 


alcohols. For water—acetone G 
vestigated G cannot be estimated owing to lack 
of equilibrium data and information about the 
degree of superheating of organic liquids. 
The value of G, however, will certainly depend 
the 


heating surface }12, 18 


ipon the size of bubbles which leave the 


and it may therefore be 


J 


exper ted to depend upon the degree of wetting 
ind the strength of adhesion of the liquid at the 


surface. If the organic component is enriched in 


1 gas-liquid interface, and also in a solid-liquid 


nterface im mixtures containing water as one 


has been suggested by some 


80). 


component 

authors [1, 2, 10 
bubbles should be expected with different types 
of organic components. However, there are not 


at present, suflicient data available to permit a 


further discussion of this point. 


a different original size of 


binary liquid mixtures 


It is found that all the maxima of maximum 
heat flux in nucleate boiling that are predicted 
by the theory are actually observed, although 
the 
‘entration of water is much less pronounced than 
the first That the 
with l-octanol than with 1-pentanol, can be 


second maximum occurring at a low con- 


one. first maximum is lower 
explained by the fact that in the former case 
the difference in boiling point between water and 
azeotrope 1s too small (see Table 1) to compensate 
the effect of superheating 


The conclusion seems to be justified that the 
position of maxima in maximum heat flux in 
binary systems can in general be predicted on 
account of equilibrium data only. 
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NOTATION 


Area of surface of test wire 


Vaporized molar fraction of liquid mixture 


Coefficient of heat transfer 


Rex 
Fig. 1) 


cal sec! em-* °C! 


Maximum nucleate boiling coeflicient of heat transfer (at C in 


cal sec"! em~* °C! 


A Equilibrium constant of most volatile component in liquid mix- 


ture 


] } and p, 


PsP. 


Heat transfer rate 


Heat flux 


Maximum nucleate boiling heat flux (at ¢ 
Solubility of component B in component A 


Pressure, critical pressure and reduced pressure 


atmosphere 
cal sec”! 
cal sec"! em-* 


in Fig. 1) cal sec! cm-* 


g/100g, weight % andmole % 


Temperature of surface of the test wire ( 


Boiling point of liquid at atmospheric pressure Cc 


Boiling point of component 


A, B and azeotrope ( 


Temperature difference between surface of wire and boiling liquid 


Temperature difference between dew point of vapour bubbles and 


boiling point of original liquid 


Molar fraction of most volatile component in liquid 


Original molar fraction of most volatile component in liquid 


(before vaporization takes place) 


Molar fraction of most volatile component in vapour 
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Abstract 
optimum temperatures profile in a 
contact times are developed 
systems. The problem is shown to 
equations 
methods 
paper numerical calculations obtained with 
The problem of two consecutive reactions A 
detail. The 


is shown in general that 


order. is discussed 
sketched It 


optimum temperature 


second in rr 


distribution is used 


Résumé On a ce veloppe 


pour un reacteur tubulaire 


tubular reactor to 
ind applications are made to 


be reducible 


and will allow the specification of the temperatures in the reactor 
an 


> B+ 


th 


appreciable gains in the 


ci-dessous les méthodes matheéeniati jues qul perni ttent 


le profil de température 


de contact minimum pour une | roduction donne 
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On montre que probleme se reduit 
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non-linéaires. (Ce 
du 
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et les méthodes a e 


la 


ten per ature 


des 
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Une publication future donnera 


analogique (REAC) 

l. INTRODUCTION 
Successful reactor design is de pe ndent upon the 
of 


pressure, 


proper choice several operating variables, 


concentra- 
be 


equivalent, the proportion of inerts in the feed, 


namely temperature, and 


tions. Pressure, or what may sometimes 


may be varied to change reaction rates or dis- 
An « xample of the proper choice 
s afforded 
product 


place equilibria. 


of concentrations | in those processes 


in which an unwanted given by a 


reversible reaction is recycled. If its concentration 
is maintained at the equilibrium value, no new 
waste product will be formed. 


* Present address 


a la résoluti 
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it 


ver pour resoudre 


ti 


des exemples numériques résolus 
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In this paper the mathematical techniques necessary for the determination of the 


insure maximum yields or minimum 


reversible and consecutive reaction 


to a system of ordinary non-linear differential 


The solution of these differential equations can be made by conventional numerical 


In a succeeding 
computer (REAC) will be presented 
the of first 
on more complicated problems is 


if the 


inalogue 


in which reactions 
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cle 


donnant la production maximun 


determiner 
ou le te mips 


Ces meéethodes sont ippliquees a des systémes 
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Because of the ol 
to temperature, the control of temperature may 
be For a 
reaction the optimum temperature will usually 


sensitivity many reactions 


very important. single irreversible 
be chosen as the highest temperature compatible 
with the materials of construction, the properties 
of the « atalyst, the problem of preheating to high 
temperatures without decomposition, and the 
problems of heat transfer which arise at high 


In 


the reaction is surface cataly sed, an increase in 


reaction rates. some cases, however, when 


the temperature above a certain limit precipitates 


- 
gas 


desorption, and this factor may well offset 
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Thus 


output value 


between installatio ists and 
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solution to both of these problems is intrinsi ally 
the to the 


profile. 


same, and leads sume temp rature 


The S¢ problems will he discussed im this ind 


It i that h 


to 


the following paper. Ss apparent Sut 
lead 


analyses. The equations which arises were solved 


the REAC, 
manufactured by the Reeves ¢ ompany 
to illustrate 


probk ms will comple x mathematical 
computer 
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the 
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on an electronic analogus 


Since 


purpose of this paper is, largely, 


mathematical methods and techniques 


physi il model assumed will he i very simple 
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1m luded sO that the 
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whole is not meaningless 


tory one in which the mathematical methods and 


numerical techniques are presented. 


* 


” « 


~ 


‘ 
‘ 


REACTIO 


MPLI 


wh h 


the 


Simpl reactions are those in 


rate of the 


the 


reversible reaction determines 


Examples of these are ammonia 


pro ess 
synthesis and sulphur dioxide oxidation 
the 


Suppose 


re action 


is reversible, second order in each direction, 
exothermic, and is carried out in a tubular reactor 
In this case the equilibrium constant decreases as 
the temperature 
high vield at the 
On the 


inlet equilibrium considerations are not important, 


increases, and it is clear that a 


reactor outlet demands a low 


temperature. other hand, at the reactor 


and the temperature may be very high in order 
rate of the reaction. Thus on this 
it that 


perature gradient is necessary to produce a high 


to increase the 


simple example appears a falling tem 


yield in a short time. It is obvious, of course, 
that 


temperature should be chosen at a low level in 


if the total contact time is immaterial the 


order to give a high yield. 


* The are indebted to K. G. Denpicu, who 


introduced them to these problems while he was a Visiting 


authors 


Professor in the Chemical Engineering Department at the 
University of Minnesota in 1953-54. 
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If G is the total mass flowing through the 
rea€tor, z is the reactor length, C, the con- 
centration of A in moles per unit mass of reaction 
mixture, a the initial value of C4, and p the 
total pressure, M the molecular weight, and T the 
temperature, then the kinetic expression for the 
reactor 1s 


M\? ' 
= (FE) [eC — kaa — C4") (0) 


> oe 
Integration at constant T and p gives, K = ra 
‘2 
M\? 2 KT? Cc K —1)+ 
fk r(* - ’ log ANY . ) tes (2) 
R/ G ky *"\C4(VK+1)—ea 
In order best to illustrate how the function F 
depends upon the temperature a simple numerical 
considered. If one chooses the 


case will be 


numerical values 


0-80 a at the outlet 
a = Cz, at the inlet, c =O =d 
1 atm. 


Mol. wet. A 54, M 


P 4-85 


46, Mol.wgt. B 
10"*, E, = 25,000 
10"*, E, = 380,000 


mn AH =5,000 
RT 


Po = 7-42 


K = 5-86 x 10% exp | + 
then the function F may be plotted as a function 
of the temperature. This has been done in Fig. 1, 
Thus, to 
design a reactor for minimum volume, one must 
operate at 550°K, and note that the minimum is a 
sharp one. It is also interesting that the optimum 


and a minimum in the curve results. 


isothermal temperature, optimum in the sense of 

fixed conversion in minimum time, occurs within 

5°K of the outlet equilibrium temperature. 
Equation (1) may be written in two forms 


first may be integrated to give 


CA 
-T*dC, 


2ax (pM\* aa. a 
( - k, C4 —k,(a—C,)* * 


G\ R; 
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Determination of the optimum isothermal tem- 
A+BVC+D, 


exothermal and second order in each direction. 


Fic. 1. 
perature for a reaction of the type 


In equation (5) the value of F, and hence z, 
depends upon how the temperature variation 
along the reactor is chosen. If one specifies how 
T depends upon z or C, then from equations (3) 
and (4) one may compute how C, depends upon 
x, and the integral in equation (5) may be evalu- 
ated. The value of this integral thus depends on 
the temperature variation along the reactor, and 
it is interesting to ask what temperature variation 
the 
This is an elementary problem in the 


will make value of this integral the least 
possible. 
calculus of variations, and Euler’s equation is the 
that 


extremal, which is in this case 


necessary condition the integral have an 


> k CY (a Cc, k. 


In section 6 it will be shown that the problem of 
calculating the maximum yield in a fixed time 
has the same answer, i.e. the optimum gradients 


necessary in each case are the same. 
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Equation (7) determines the temperature if the 
concentration is known and hence the temperature 
is fixed locally and does not depend upon the 
history of the system. Equation (6) says, in 
effect, that the reaction rate is to be maximized 
with respect to the temperature at each point. 
The temperature variation along the reactor may 
then be computed using equation (7) and a step 
by step procedure using equation (8). From 
equation (7) at the entrance of the reactor, 
C, =a, T = E,/2R. This temperature is much 
too high in general for practical operation, and 
some arbitrarily fixed temperature is 
chosen. If one chooses, for example, this upper 


temperature at 870°K, operating iscthermally at 


upper 


this temperature until the concentration and 
870°K satisfy equation (7), and follows this 
section of the reactor with a temperature gradient 
calculated from equation (7), it is found that 
F = 805. For the isothermal case F = 648, and 
therefore the reactor volume if the 
optimum gradient is used. It should be pointed 
out, however, that heat removal requirements 
have not been considered, and could be so severe 
as to make the optimum gradient impossible of 
attainment. 

Problems similar in nature to this one have 
been treated by Densicu [4], ANNABLE [1], 
CALDERBANK [8], and Temxin and Pyzuev [7]. 


is halved 


8. THEORY FOR TWO 
CONSECUTIVE REACTIONS 


Consider the system of two consecutive reactions 
A+B--+C 


The feed is component A, and B is the product 
of interest, the yield of which is to be maximized. 
Component C is a waste product. 

If one supposes that the activation energy of 
the second reaction is greater than that of the 
first, then a high temperature will produce more 
C than desired. It is therefore apparent that, 
near the end of the reaction time when the con- 
centration of B is high, the temperature must be 
low. On the other hand, at the beginning of the 
reaction when the concentration of B is low the 
temperature may be high, since there is little B 
in the system and hence little C can be formed. It 


seems intuitive therefore that if E, > E, there 
should be a falling temperature gradient along 
the axis of the reactor. 

The method of attack will be demonstrated on 
the simple case in which both reactions are of the 
first order. The technique will then be extended 
to a more general case, more interesting prac- 
tically. The letters z and y in the future will be 
concentrations and ¢ will be the time. Since 
pressure changes are to be neglected, the equations 
for the tubular reactor have essentially the same 
form as those for the batch reactor. Assume then 
that the kinetics of the system are described by 


a= —a(T)z (8) 


Y = a(T)2—b(T)y (9) 


where a(T) and 6(7) are the specific rate con- 
stants. Our aim is to determine the temperature 
T (t) and therefore the temperature profile in the 
reactor which will give the maximum value to 
the yield y at a given reaction time ¢ equal to @. In 
the course of the problem a(t) and 6 (t) will be 
the specific rate constants obtained by replacing 
T by the appropriate function of the time. The 
solution to the set of equations (8) and (9) is 
given by 


t 
-fe (B) 48 
0 


z(t)=a,¢ 


‘ 
—f > ae 


y(t)=yoe° 


‘ Seinen fener 
+2,{ a(p)e e” 


dp (11) 


where z, and y, are the feed concentrations of A 
and B. The expressions z(t) and y(t) describe 
the composition of that part of the reacting 
mixture which entered the reactor at time zero. 
At every instant ¢ during the course of the 
process, the part of the reacting mixture under 
observation assumes a temperature 7. Thus let 
T, be the temperature at the instant ¢,, with 
0<t,<@. The yield of y thus depends on an 
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; at the 
The 


function y may then be considered as a function 


infinity of variables, the temperature T 
instants ¢; during the course of the process. 


of the independent variables T;. If one supposed 


4° 
that the reactor were divided up into a finite 
number of sections, each section having a tem- 
perature 7',, then equation (11) could be written 
as a sum of terms, each term describing a certain 
isothermal section of the reactor. The yield y at 
the reactor outlet would then be a function of all 
these temperatures and the condition for an 


extremal would be 


dy (8, T;) 
dT; 


0 forall: 


However, here the variables 7; are not discrete 
but the 


interval! (0, 0). The dependence of y on T is then 


form an infinite set defined on time 
described mathematically by saying that y is 
a functional of the temperature, a functional being 
a quantity which depends on the entire course of 
one or more functions rather than on a number 


of discrete variables. The notation 


y (9, [T}) 


will be used to denote this functional dependence 
of the yield on the temperature. 
Consider a function 7, (¢) which gives rise to a 


y (0, [7,(t)]) and a function 7, (t) 


which is equal to 7, (t) everywhere outside an 


yield wo = 


In this interval 


T, (t) + 8T (t) 


interval t, <t < t,. 
T, (t) = 


The yield with the temperature 7, (t) is y, = 
y (6, [7 (t)]) and the average variation in y due 
to the disturbance 7' (t) around t = ¢; is given by 


% 


Ay (@,{T]) = [y(9,[7(t)]) — y(@,[To(t)] ) ae 


ta 
ir) % 


1 fo, mnaso(ferurs) 
la ta 


with At = t, 
series of the functional y. 


t,, after a development in Taylor 
This may be written 


1 dy (8, (TKt 
ay (9, (7) = % Bal we 


% 
[ ST (t) dé 


+ 


where t, <t; < ty. 


P As t, >t, and t, +1t,, then 


the point variation of the functional y may be 
defined as 


Ay (6, 
By (8, (TVery, = lim —2 a) 


Calling 


(18) 


expresses the dependence of the yield on the 
temperature encountered by the reacting mixture 
time ft t 
partial functional derivatives. 


at a Such expressions are called 


” 
They have been 
introduced mostly by VoiTerra [8] and a good 
account of functional analysis may be found in 
Levy [5, 6]. 
be expressed by setting the partial functional 


Hence the extremal condition will 


derivative of the yield with respect to the tem- 
perature equal to zero for all values of ¢; in 
(0, 8). This condition, which is essentially equa- 
(12) will then 
determine the values of the variables T',, that is 


tion extended for functionals, 
to say the function T(t). which describes the 
optimum temperature profile in the reactor. 

In calculating the quantity given in equation 
(18) it is simpler, intuitively, if one thinks of the 
T's as being discrete variables. In equation (11), 
B varies from 0 to @ and ultimately passes the 
point ¢, at which the variable 7; enters into play. 
The partial functional derivative will therefore 
have different contributions arising from the 
different portions of the time interval (0, @), 
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Carrying out this operation, there results 


a 
—f o(7) ap 


MOLT) _ _ gor) ye ® 


oT; 


8 a 
—fa(T)4y —fo(rdy 
8 dp 


uy 
—bxT to { a( Tye 
0 


. 
(T)dy 


i 
-f a(T)dy - f 


e 0 
+ az(T;) ae . 


—ar(T,)to [ (Te ’ 


* 


4 


. 8 
-f a(T)dy - f b(T)dy 
B 


dp 


The extremal condition, obtained by setting 
this equal to zero for all i, becomes 


a7(T;)2(t,) = bx T,y(t,) — 


B 4 
—{ a(T)4y f o(r)4y 
F) 


To A7/ r)| a(T je ' 


i 


dp 


As this expression must be valid for al! values of 
t, between 0 and @, one may replace ¢, by t. 
Replacing also T by its unknown expression as 
a function of the time, and setting 


(18) 


one finally obtains P 
-f a(y)dy —[ %(r)¢y 


‘ "dp (14) 


- 
a(t) — A(t)y(t) <1 a( Bye 


In order to obtain a more manageable expression 
equation (14) may be differentiated with respect 
to the time to give 


Ap(T) = [ocr) —a(T) A (n)| = (15) 


Equation (15), when combined with a suitable 
auxiliary condition, is completely equivalent to 
equation (14). 

Thus the solution of the system of simultaneous 
differential equations (8), (9) and (15), gives the 


86 


solutions z(t), y(t), and T (t) dependent on the 
initial values 2, yp and 7, which are related to 
@ by equation (14) and which may be put in the 
form 


6 
fe (2. Ve» Te ¥) 


0 

z,+[{1—A(T,))y, = We dy (16) 

If z, and y, are given, equation (16) determines 

the relation between the process time @ and the 

optimum initial temperature 7. Formally then 
our problem is solved through the system 


—a(T)z (8) 


dt 


dy 


T —)5(T 
77 a(1)z (1) y 


(9) 
dT _b(T) —a(T) A(T) 


ao AT) 


qi (15) 


The kinetic constants a(T) and 6(T7) are 
usually of the Arrhenius form and the above 
equations are non-linear but may be solved 
numerically by well-known methods. Such a 
numerical solution may be performed if the 
initial values z,, y,, and 7, are given. However, 
the quantities usually given will be z,, y, and @. 
In this case one must assume an initial value for 
the temperature, 7,. The differential equations 
may then be solved. The time for which this 
solution satisfies the extremum conditions is 
given by equation (14), for when t = @, 


z(@) — A(@)y(@) = 0. 
For convenience 


Z(t) =2(t) —A(t) y(t) (17) 


or similar functions which will be encountered 
later, will be called the “ zero function.” Z (t) has 
the property of being equal to zero at t = @, 
and @ is thus easily determined by following the 
behaviour of Z (t) during the course of a numerical 
solution starting with an assumed value of T>. 
The temperature profile itself will be obtained by 
eliminating z and y from equations (8), (9), and 
(15) to give 
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aT ; a(T) — B;(T) —4y 
df? B(T) dt 


[a (T) b (7) =o (18) 


where 


b T) A(T 
B(T) (T') nee (T) (18’) 
=F" 


The initia] conditions to be associated with 


dT 
equation (18) are 7, and » and are related 
dt 


0 
to %, Yo, and @ through equation (16) and 


equation (15) 


To 


iT 
¢ == B ( T.) 
0 : 


dt 


40 


This last relation shows that if there is no 
trace of component B in the feed, the optimum 
solution calls for an infinite temperature gradient 
at the reactor entrance. It should also be remarked 
that rather large deviations from the theoretical 
optimum temperature profile should have a small 
effect on the optimum yield of B because of the 
This 


effect is counterbalanced by the sensitivity of the 


stationary nature of the optimum profile. 


rate constants to temperature. 


4. A MoRE CASE 


To illustrate further the method of derivation of 
the equation defining the optimum temperature 


COMPLICATED 


profile, a more complicated case involving both 
consecutive and side reactions will be treated. 


Consider the kinetic scheme 


>B+xX 


‘ 


> Y 


Feed component A gives the product of interest 
B by a second order reaction. Product B decays 
to waste X by a first order process. At the same 
time a competitive reaction occurs ; feed com- 
ponent A is destroyed by a first order process 
It is desired to 


this model the 


and gives a waste product Y. 
maximize the yield of B. In 
kinetic scheme is 


dy 


(20) 
dt 


=a(T)2* —b(T)y 
where a, b, c are the specific rate constants of the 


> B, B > | & A- ar or 


modifications if volume changes 


respective reactions A 
appropriate 
during the course of the reaction are serious. 
Equations (19) and (20) may be solved in closed 


form to give 
(21) 


-[d(r)dy 
. 0 
TX x,a(8)H(B)*e (22) 
0 
where H (t) is defined by 


i 


Hity' =1 + 9 | a(B)e 

o 
As in the preceding case the temperature has been 
replaced by its expression as a function of reaction 
time in the kinetic rate constants. 

The partial functional derivative of the yield 
of B with respect to the temperature at a given 
time ¢ will be evaluated. The process is carried 
out as before, the intervals (0, t) and (t, @) giving 
different contributions to the derivative. How- 
ever, the different functionals of temperature 
which enter into the expression of the yield are 
more complex and the example provides a good 
illustration of the application of the method to 


more difficult cases. Thus 


dy (9, [7] 
 — y 
—f e(e)de 
) 


2 { tae) H(B"| ext | atre dy 


* 


6 ' 
“ 
t 


t 
—f e(o)do 


arit) e 


b(t) [a a(B) H(B)* e 


0 
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° 
-f d(y)dy 


brit) ye , 


v 
—[ indy 


+ 23 ay (t) H(t)* 


After substitutions from equations (21) and (22) 


and elimination of the common factor 


f (rar 


f 
e 


the cumbersome above give the 


following equation, which expresses the vanishing 


expressions 


of the partial functional derivative of y, 


l I 


I To 


tf 
P serra 
F ss.” . 
2a (8) ai 8)’ e dp 


* 


ar(t) 2*(t) brit) y er(t) | 


arit) (28) 


on (14) 


nperature 


Equation (23) is the analogue of equat 
The differential for the 


gradient may be obtained as before 


equation 


cAdy + 2a r* | 


CT { T) 
ar\| T) 
b,(T) 


A(T) a(,T) 


It would be the 


centrations z and y of components A and B 


possible to eliminate con- 


among equations (19), (20), and (24) in analogy 


with what was done in the previous example, but 


the complexity of the expressions thus obtained 
makes this procedure lose its interest in the 
present case. 

A remark in taking the partial functional 
If one thinks 
of the integrals appearing in the expression for 


the yield of B as sums, infinite to be sure, then 


derivative may be worth making. 


the derivative may be taken in the standard way, 
and affects only those terms of the sum which 
Once the 
derivative has been taken the sums then revert 


depend on the value i of the index. 


to integrals. This is certainly not rigorous but it 


is a convenient way of thinking about the 


problem. 


5. Tue Imeprticitr MetTuop 
The method explained above is dependent upon 
solutions of the kinetic 


It is a fact, however, 


the existence of closed 

equations of the system. 
that even in rather simple cases — for instance a 
first order reaction followed by a second order 
one — these closed solutions are not known or, if 
known, are unwieldy. This does not prevent the 
solution of the 
methods It 


devise a method of obtaining the equation which 


kinetic equations by numerical 


seemed interesting, therefore, to 
defines the optimum temperature profile, which 
would not make use of closed solutions. Consider 


again a system of two consecutive reactions 
A+»B--—+C 


where B is the 
that the 


more general form 


omponent of interest and assume 


kinetics are now represented by the 


dx 
dt, 


a(T) F(z) 


dy 
dt 


na(T) F (zr) b(T) G(y) (26) 
Thus, 


of the form 


this representation may include reactions 
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The letters z and y designate the concentrations 
of components A and B. a(T) and b(T) are 
the kinetic rate constants of the two consecutive 
reactions, and F(z) and G(y) depend on the 
order and assume mechanism of the reactions but 
The 
coefficient n has been introduced to include cases 
of the form A + 2B 
in stoichiometric 


do not contain the temperature variable. 


» C where there is a change 
like the 
second example above, where a parallel reaction 


coefficients or cases 
of essentially the same mechanism competes with 
the reaction giving component B. 
Consider then the partial functional deriva- 
tions 
dy (t, [T 
and “Y L !) 
oT; 
the index i on T being dropped whenever there is 
no danger of confusion. The 7’;’s are independent 
variables and are not functions of the time t. Thus 
the order of the differentiations with respect to 
tT; 


derivative thus obtained. 


and ¢ may be changed without affecting the 
Therefore, one may 
write, using equations (25) and (26) 


8 (t — t,) [nay (T) F (2) 
by (T) G (y)) 


+ na(T) F’ (2) Tr, 
The symbol 8 (t — t,) is the Dirac delta function, 
and no rigorous mathematical significance should 
be attached to its use. It should rather be under- 
stood as a convenient symbolism, the proper use 
of which allows us to bypass a rigorous but 
cumbersome derivation of the result. 8(¢ — ¢,) is 
equal to zero for ¢t + ¢,, and is infinite at ¢ = ¢, 
such that 


ire 


d(t —t)dt=1 


* 
'i-« 


Taking into consideration the fact that for 


0 and 


Equations (27) and (28) may be solved. Equation 
27) is an ordinary differential equation in 


whose solution is 
F 
foinP tear 
t;) ay (B) F (2) e* dg 


and hence (since 6(8 —t,;) is a kernel function 


for a singular integra!) 


i 
- fa y) FP’ (2idy 


t 
ar (T;) F (a,)e 


where 7'; and z; are the values of the temperature 
and the concentration of A at time t;. Equation 


(28) may then be written in the form 


\ 


b(T) G’ (y) 


dy 
oT; 


8 (t — t,) [nay (T) F (z) 


i 
—f a y) F’ (z)dy 
oY 


—na(T) F’ (2) a7(T) F (ae 


— by (T) G (y)} 


Then 


t 


oY | 8 (8 —t,) (naz (T) F (2) 


7; J 
0 P 
Joineundy 


~ bp (T) G(y)je dp 


— | na(T) F’ (2) ay (T;) 
-] B 
—{ aly) Playdy + feinawedy 
i 


F (ae . dp 
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i 
fo y)G'(y)dy 


— bp (T,) G (y,)] e 


i 
fone (y)dy 


—nar(T,) F (ae ; 


> ' 
- = [naz (T;) F (z,) 


B 


‘ feomew a( y) F’ (2) \dy 


a(T) F’ (z)e" dg (29) 


u 


The extremal condition on the yield of B at 


, , d ; 
time 8 demands the vanishing of y Cancelling 
) 


. ' 

the common factor, 

ij 

fone (y)dy 

‘ 

e 
replacing ¢ by @, and suppressing the index i, 
there results 


na7(T) F(z) — b7(T) Gly) - 
8 
fnew a(T) F’(2)\d 


a(T’) F'(z) e 


na7(T) F(z). 


C) 


dp (80) 


It should be remarked that T and z or y which 
enter as arguments of the functions under the 
integral signs, should be replaced by their ex- 
pressions as functions of the respective summation 
variables 8 and y. Using the definition of A (T), 
equation (13’), 


5 
pene war? (z)}dy 


[ an) F'(z)e 


dp (81) 


As before this integral expression may be trans- 
formed to a differential equation to give 

+ [a7 SY 
n dt 


a = aT) F(z) 


— Gy) 
+ [1-4 ee 


| [oe G'(y) — a(T) F'e)| 
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which becomes 


aT nB(T) 
dt 


F (x) G'(y) (32) 
G(y) 

where B(T) is defined in equation (18’). An 
immediate check is provided by a comparison 
of equations (15) and (32). 

dT @T dT 
dt’ de® de® 
would give the differential equation defining the 
shape of the temperature profile. Such 
elimination may be carried out if F (2) and G (y) 


An elimination of z and y between 


an 


are given in explicit form, but the results are 
usually too complicated to justify the trouble. 


6. EQUIVALENCE OF Two PROBLEMS 

In the introduction it was stated that another 
approach to the study of the optimum tempera- 
ture profiles was possible ; namely, to search for 
the minimum process time @ at which a given 
yield y,, of the product of interest can be attained. 
This problem is difficult to treat by the direct 
approach, but falls well within the possibilities 
of the implicit method. 


application of this method, the problem for a 


As an example of the 


system of two consecutive first order reactions 
will be solved. Consider again the system 
A + B +C for which the kinetics are given by 
equations (8) and (9). The yield y of the com- 
ponent of interest B will be taken as the inde- 
pendent variable, and the reaction time ¢ and 
concentration of feed component z as functions 
t(y) and z(y). This representation is proper and 
will not to multivalued functions, 


give rise 


provided values of yup to the maximum (3 = 0) 
and not beyond it are taken. | 
Equations (8) and (9) may be written 
—_ 1 — = 
 a(T)a —b(T)y 
dz —a(T)x 
dy a(T)x—6(T)y 
The process time @ is given by 


A 
z—b(T)y 


J am 
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The limits y and y,, are fixed, and T and z are 
to be replaced by their expressions as a function 
of y. It is desired to make the process time @ an 
extremum with respect to the variable 7',, which 
is now the temperature at y= y,. The partial 
functional derivative of 8 with respect to 7, is 


| 


To calculate the partial] functional derivative of z 
with respect to 7;, apply the implicit method to 
equation (84) 


l 
(az — by)* 


or 
a 
dT; 


(azz - brw)| dy (86) 


+ 3(y —y,)- 


dw ~~ aby Mw . 
dy2T; (ax—by)* dT; Ky—y 


This gives 


ow 


be, —ab, 
iT; . 


(az—by)? 


0 


| ys 
ry e 
vow 


The summation variables 8 and y under the 
integral signs stand for the yield y. If equation 
(87) is substituted into equation (36) and the 
delta function term is integrated, one finds 


} (ax — bBy* 
v 
The extremal condition on @ is therefore 
art — bry 4 
(bay — aby) zy 
vm 


J 


a 


(az — bay?” 0 


(38) 
On taking the corresponding differential equation 


load] 


dy 


a 
(az — by * 
2 
(b —aA)ay (az — by)? — 


_a— Ay 
(6 —aA) zy 


0 
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Since 


there results 


which is the same expression as the one which 
was obtained for the temperature gradient for the 
problem of the maximum yield at constant pro- 
cess time. It may be remarked also from equation 
(88) that the zero functions are the same. Thus, 
at constant process time @, there is a maximum 
yield y If one 
seeks the shortest process time ¢ in which to 
- 8. The two 
problems are therefore mathematically equiva- 
lent. This has been proved for this one case 
only, and to be rigorous should be proved for 


y,, of component B obtainable. 


obtain a yield y,, one finds that ¢ 


more gencral cases. However, from the calculus 
of variations it is known that if an extremum 
of a functional 1[T] under the condition of 
restraint 


constant 


J[T) 


is sought, it is necessary to find the extremal 
functions which satisfy 


bl + pid = 0 


according to the method of Lagrange multipliers. 
But this would also be the condition which would 
have to be satisfied if the extremum of J [7] is 
required under the condition of restraint 


I [T] = constant. 


Therefore the result found might have been 
expected. More details about this may be found 
in Courant and Hi.Bert [2]. 
7. DISCUSSIQN 

In this paper methods of determining the optimum 
temperature gradient in a tubular reactor have 
been developed. The problem has been reduced 
to the solution of a system of ordinary differential 
equations which will be discussed in a succeeding 
paper. No rigorous justification of the procedures 
used has been given. While this seemed jutsified 
during the analysis in order to maintain its 





Ovecu Britovus and Neat R. AMUNDSON 


practical value, some discussion of the points 
under question seems proper at this stage. 
The first question which naturally arises is the 


validity of the solution. Will the temperature 


profile as calculated provide the maximum yield, 


or is there another arrangement of temperatures, 
for instance operation at the highest attainable 
temperature, which will give an even higher 
The question also remains unanswered 
whether the problem has an optimum solution at 


yield ? 


all. This sort of question always arises when the 
calculus of variations is used. To be complete 
the that the extremum 
obtained is and further that it 
is the maximum maximum in the sense that there 


proof should show 


a& maximum, 


is no other temperature distribution which will 
give a higher yield. The first part concerning the 
existence of the maximum can presumably be 
answered by examining the sign of the second 
variation while the second part requires variation 
in the large, rather than the small. These demon- 
strations are always very difficult, and it iscustom- 
ary in physical problems not to carry them 
through but to depend on intuition and numerical 
This 


is unsatisfactory to many, but should be accept- 


procedures to rule out the doubtful cases. 


able to most engineers and applied mathema- 


ticlans. 


NOTATION 


A, B,C, D, ete. 
b’7 (T) 
a’7(T) 
a(T), b(T), c(T), ete. 
stants, generally of Arrhenius form 
b(T) — a(T) A(T) 
. A’ (T) 


= chemical specie taking part in a reaction 


specific reaction velocity con- 


eT ( T) 
— ay (T) 
= activation energy of a reaction 
= function defined by equation (2) 
= concentration effect in kinetic expression 
= concentration effect in kinetic expression 
= specific velocity constants 
= equilibrium constant for reaction 
= 29 + Yo + %, sum of initial concentrations 
gas constant 
activation energy ratio in a reaction 
temperature in reactor 
temperature at time /; after entrance 
time 
particular instant of time 
concentrations of A, B, C, respectively 
zero function—differently defined in each problem 


' 
z,y,2 
Z(t) 


@ = minimum time 


Prime denotes differentiation with respect to subscript - 
generally T. 
Subscript zero denotes initial condition. 


Subscript 1, 2, 3, . . . denotes a particular reaction. 
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The application of linear programming to design problems 
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Abstract 
profession a recently developed mathematical technique, called “ linear programming”’, 
is exceedingly valuable in solving certain general problems connected with, inter alia, the efficient 
design of a chemical plant, and the discovery of an optimum production schedule which would 
The general and more important features of linear pro- 


It is the purpose of this paper to bring to the attention of the chemical engineering 
which 


maximize the profits of an enterprise. 
gramming will be discussed in the main body of this paper, while, in the Appendix, an illustrative 
problem of, it is hoped, general interest will be worked out numerically. 

Résumé—Cet article a été écrit pour attirer l'attention des Ingénieurs du Génie Chimique 
sur une technique mathématique récemment développée et appelée “ linear programming,” 
extrémement précieuse pour résoudre certains problémes généraux relatifs, entre autres, a 
lobtention d'un projet efficient d’une installation chimique et & la découverte d'un programme de 
production maximum qui porteraient & son maximum les profits d'une entreprise. Dans le 
corps de cet article, les auteurs traitent des caractéristiques générales les plus importantes du 
“ linear programming,” tandis qu’en appendice ils étudient numériquement un exemple d’intérét 


général 


A SIMPLE PROBLEM 


; = rate (or level) at which process i is oper- 
1, 2). 


amount of raw material 7 consumed in 


It is believed that an easy introduction into the ated (i 


main ideas of linear programming can be made by 
process i, when the process i is operated 


considering the following simple problem. 

A certain entrepreneur has at his disposal two 
raw materials, A and B, fixed price, which can be 
reached in different ways to produce two propucts, 
C and D, which he can sell, again, at a fixed price. 
We shall assume, moreover, that the two pro- 
cesses, which the entrepreneur may employ, are 
fixed in the sense that both the amounts of A 
and B consumed in, and the cost of operating, 
each process per pound of product, are pre- 
determined parameters which cannot be altered. 
The following then If we 
suppose that the available supply, per unit of 
A and B cannot exceed, for one reason 


question arises : 


time, of 
or another, certain upper limits, then at what 
rates must the entrepreneur produce the two 
products, C and D, in order to maximize his 
profits ? 

The above problem will now be solved by an 


elementary method. Let first 


1, 2). 


net profit from selling the product of 


at a unit rate (7,1 


process i, when this process is operated 
at a unit rate. 
maximum available rate of supply of the 
raw material ¢. 


r = the total profit per unit time. 


In mathematical language then, it can be said 
that the entrepreneur desires to maximize the 


linear function 

r= C2, + gas (1) 
where the only two variables at his disposal, 2, 
and @, are subject to the conditions 
(2a) 
(2b) 


Gy, 7; + Ayg% = & 


Ag, 7 + Age Te > 89 


and where the parameters a,,, 8, and ¢, are fixed. 
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Morcover. we must obviousls re that, 1. T pound of uf | pound of B + 2 pounds of ¢ 
2. 1 pound of 2 pounds of B + 3 pounds of D 
ry 0 and wv, 0 (2c) 


and 
Now. let us consider the quadrangle OGHK, a 


100 pounds of 4 day 
shown in Fig. 1 It is at once clear that any 


point (27,.2'.) which satisfies equations (2a), (2b) 150 pounds of B day 
must not lie outside OGHK, and more- 1 dollar per pound of ¢ 
ver that the line me 0.80 dollars per pound of D 


Therefore equations (2a) and (2b) become 


] 
Ze. < 100 
3 2 


150 


’ ‘ 4 . 
Now, since ; » the point H, in Fig. 1, 
Cy : 


represents the optimum production schedule, and 
z, = 100 pounds per day 
7, = 150 pounds per day 
On the other hand, if 
C, = 0-60 dollars per pound of D, 
then the optimum production schedule would be 
\« z, = 200 pounds per day 





ty = 0 


Fic. 1 ; 
We perceive therefore that, in this simple case 


pn in any event, the production schedule, which 
must pass through OGHK, or one of its four : 


; maximizes the total revenue r, depends in a 
corners. Finally, since r is the only unknown 


discontinuous manner on the various fixed para- 


parameter in equation (3), the maxiumm value 
meters of the problem, namely the coefficients 


of r is proportional to the greatest permissible 

wh-ge-guean . — . a,,, ¢, and s,, and that, moreover, this optimum 
intercept of equation (3) with the z,-axis. It is 

, point (z,,2,) is situated at one of the corners of 

the convex polygon OGHK. Finally, it can be 


said that, due to the peculiar nature of the 


not diflicult therefore to show that three cases 


the point (z,,.7,) which restrictions which the two variables z, and 2, 


must satisfy, this problem just illustrated could 


IMAXUINIZCS TF, 
not have been solved by the methods of differ- 
~ Aya ‘ : > 

2. if = « then the optimum pro- ential calculus. 
Qi Cy ag, 


cdiucttw *( s,s 7 *! by 9 ° . 
ction schedule is given by H Tue MATHEMATICAL FORMULATION 


a OF THE GENERAL PROBLE™M 


12 


Cc. 
Ss. If 3. 
‘ Gi) 


» then G is the optimum point. 
So far, only a very simple special case has been 

{s a simple numerical example consider the considered. However, the problem can easily be 
following: The two reactions are generalized as follows. It is supposed now that 
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the 


materials, in limited supply, which he can use 


entrepreneur has at his disposal n raw 


ink pre defined processes to produce certain goods. 
And again he desires to determine the optimum 
production schedule which would maximize his 
profits. Stated mathematically, one is asked to 
calculate the values of z,, 7,....2, (that is, the 
rates at which the various processes should be 
which would the linear 


operated) maximize 


expression 


(4) 


To 7 


C12, + Ce%2 


with the restrictions 


G5, 7, T Gjg%2 


and 


the 
assumed to have been specified a priori. 


Of course, parameters ¢;, Gj, 


‘ and 8; 


Now, it was shown earlier how a simplified 
version of the above problem could be solved by 
a graphical method. More complicated cases can 
be handled by a recently discovered technique 
called 


linear programming has been defined as the study 


linear programming. Strictly speaking, 
of the maximization or minimization of a linear 
expression, where the variables in that expression 
are restricted by a set of linear inequalities. Such 
a system is described by equations (4), (5a) and 
(5b); it can be solved by a number of methods, 
the most common of which is the “ Simplex 
Method.” 


an iterative procedure ; 


In essence, 


this “Simplex Method ”’ is 


that is, one starts with 


a set of values for 7,,2,,....2, which satisfy 


(5a) and (5b), and then calculates 


of values for the above mentioned 


equations 
another sct 
variables, which increases the total revenue. This 
procedure is then repeated until r is maximized. 
In addition, the method is simple enough, to be 
adaptable to machine calculations. 

Problems in maximization, or minimization, 
which can be treated by the methods of linear 
programming, are encountered very frequently 
in connection with, to mention only a few cases, 
the efficient design of a chemical plant, and a 


large variety of planning operations. To be sure, 


For 
one thing, the assumption that the purchasing 
price of the raw materials and the selling price 


there are severe limitations on the method. 


of the products are known in advance and are 
independent of the rate of consumption or the 
rate of production of the substances is, usually, 
an over simplification of the real situation. Also, 
the *“ Simplex Method,”’ in its present form, holds 
only when the number of processes available to 
the entrepreneur, is finite. In spite of these 
restrictions, however, there are large numbers of 
to the 


system of equations (4), (5a) and (5b) would 


important practical problems, which 
apply with sufficient accuracy, and which there- 
fore could be investigated by the methods of 
linear programming. For a much more detailed 
account of the subject, as well as extensive 
references, the books by 
DorrMaNn [1] and by Carnes, Cooper and 


HENDERSON [2]. 


one could consult by 


APPENDIX 


An ILLUSTRATIVE EXAMPLE 


The numerical example to be presented now-will be worked 
out by the so-called ** Simplex Method,” which is the most 
commonly used procedure for solving problems of the 


f 


linear programming type. The notation of Reference [2] 


will be used throughout this appendix. This reference 


should also be consulted for the details of the ** Simplex 
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Fic. 2. Flow diagram for the illustrative problem. 
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Method,” as well as for instructions concerning the use 
and construction of the “ Simplex Tableau." The purpose 
of this example is merely to illustrate the types of prob- 
lems that can be analysed by this method, and to give 
some idea about the amount of work needed to perform 
the calculations. 

The flowsheet of a hypothetical chemical plant is shown 
in Fig. 2. The plant produces four products, G, BE, F and 
H, and consumes four raw materials 4, B,C and D in 
limited supply. The four available processes are depicted 
schematically in Figs. 8 to 6. 


SO 
ol 


Frio. 6. 
Processes 
4 and 5. 


Schematic diagrams of Processes 1, 2, 8, 4 and 5 of the 
Ulustrative problem. Processes 4 and 5 are identical in 
all respects except for the net profit. 


Table 1. 


The following data are given : 
Raw material Mazimum available supply 
400 Ib. /day 
800 Ib. /day 
100 Ib. /day 
250 lb. /day 


Pounds of raw 


Cost 
$1.50/Ib. 
$2.00/lb. 
$4.50/Ib. 
$2.50/Ib. 


Cost of 
reaction 
$1.50/lb. A 
$0.50/lb. A 
$1.00/lb. G $3.75 /\b. 
$2.00/lb.H $5.00 /Ib. 


It is also stated that the reaction cost of Process 4 is 
$2.00/lb. H only up to 75 lbs. of H per day, and that for 
any additional H produced, the reaction cost is $2.20 per 
pound of H. So, in essence, the number of processes is 
equal to five. 

Following a well established procedure [1], we shall 
introduce five disposal processes, numbers 6, 7, 8, 9 and 
10, which dispose of, respectively, raw materials A, B, 
Cc’, C’ and D. C”’ is that part of the substance C which 
is used exclusively in process 4, while C’’ is used exclusively 
in process 5. 

If now, 2, 1 qi <5, denotes the net profit from 
process i, in dollars per day, and z,, 6 <i < 10, measures 
the rate of disposal by process i in pounds per day, the 
following five equations are obtained : 


Selling price 

of product 
$8.00 /Ib. 
$2.884/Ib. 


8 1 2 
2, + Bg + 55 + 5%e + se + Weg + Oty + Oty + Ory +020 = 400 
lz, 4 Ly + Sty + stg + 2 + Oty + Ly + Oey + Ory + 02,0 = 800 


Oar, + Oaty + Oaty +g +0ity +Orrg + Osty + lity + Ory + Org 15 


Simplex Tableau, First Stage. 





P; P, | 





1/8 1/2 | 





-1/3 





0 

















1/2 











5/6 





























Min “2 = BO ati =m 6 = 4, ag = a om 
i ™% 


2u4 ia ty yy — yt 





The application of linear programming to design problems 


Table 2. 


Simplex Tableau, Second Stage. 
































Zz. 
Min -2° 


Now use z, as the new active process 
Lis 
3 J 


Table 3. 


Simpler Tableau, Third Stage. 


























Profit 2484 











Or, + Ostg + Ost + Oar, + Lrg +Orrg + On, + Org + ltg+Or, p= 85 


1 
ery + Oity + Leg + 5g + Lrg + ety + Ory + Ory + Oey + 12,9 = 250 


The first “Simplex Tableau” is obtained by letting 
fg 2g — 24 = 25 = 2, = 0, and the procedure there- 
after is identical with that employed by [2]. It will be 
seen that the problem will be solved by just three itera- 
tions. 

We see then that the “ Simplex Criterion ” [1] is now 

16 
satisfied, and that therefore r = 907; dollars per day is 
the maximum attainable profit under the conditions of 
the problem. Four active processes are used, namely pro- 


s 
r ; 
9 10); 


cesses 1, 3, 4 and 5, together with one disposal process, 
No. 9 
2 give the same profit, the levels at which they operate 
varying the total net 


It can also be observed that since processes 1 and 


could be interchanged without 


revenue 
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Table 4. Simplex Tableau, Fourth Stage. 
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Book reviews 


A. Hor 


Process Principles, 


en, K. M. Watson, R. A 
Part I, 2nd 
Hall 


Chemical 
04 Py 
oe \ 68s 


prin 


nis is a second editior of 
d t 


of the few books to treat specific ally the 


an important b ok that can be 


warmly recommende all chemical engineers. It is one 


subject of material 
The ject with all 
should both 


because of its utility in drawing up quantitative flow sheets 


ind energy balances sul which 


be 


is one 


chemical engineers intimately familiar, 


for plant design, and its value in interpreting experimental 


results from existing plant. These two applications are 
amply recognized by the authors who give large numbers 
of worked examples and problems from industrial practice 
The book but 


sufficient modern data on vapour pressures, heat capacities, 


is primarily a University teaching text, 
heats of formation, combustion and dilution are given to 


ike it valuable as a reference work for the design or 


plant enginecr 


The second edition has two more chapters than the 
first, one dealing with mathematical procedures and the 
The of the book has 
revised, new problems have added, 


numerical answers are now given to about half the prob 


other with adsorption rest been 


carefully been and 


lems, and ought to have been given for the lot ; numerical 
problems without answers are almost useless to the student 
working by himself. The printing and layout have been 


much improved ; unfortunately the price has increased 


The book could usefully be shortened by deleting the 


S.W. Ferris : Handbook of Hydrocarbons, Academic 


Books Ltd., London, 1955. 324 pages. 68s 


Tut writer of the book, who is affiliated with the Research 
Development Department of the Sun Oil Company, Marcus 


Hook, Pa., 


The book is indeed a help to those investigators who are 


U.S.A., has written a valuable contribution 


connected with hydrocarbon research. 


In Table A the hydrocarbons are listed in the order of 


the boiling points 
index), d (specific gravity), and melting point are included. 


at 760 mm of mercury; n (refractive 
In Table B the listing is done in the order of the empirica! 
formula. The hydrocarbons of the same type are brought 


in separate groups. 


Both tables take up a total of 200 pages. The remaining 


$4 pages are used for Table C which gives alternate names 


129 


99 


first three chay ters These deal re spec tively with mathe- 


nrc 


matical procedures, stoichiometri 


relationships, and the 
behaviour of ideal gases 
thi 


rt 


The student who has not grasped 
elementary mathematics, chemistry and physics will 
n read a bor chemical 
he 
the properties of triangular diagrams, excellently explained 
1, could be 


systems, rather than treated in the abstract. 


k on process principles, nor will 


learn it from so few pages. Specialized matters such as 


in Chapter moved to the part dealing with 
ternary 


Another addition which has a doubtful place in this 
box 


is a somewhat detailed study of stagewise solvent 
This would seem to be 
of 


kind would lead to a book on Unit Operations under 


extraction processes in Chapter 6 


trespassing on another field, and further excursions 


this 
the title * Material and Energy Balances.” 


The new Chapter 10 on adsorption is a valuable expan- 


sion of the short section in the first edition. There are 


numerous which have cost 
the 
of 


New examples of current industrial interest are a 


also revisions must much 


25°C 


Modern data on critical properties is 


thermochemical data 


18°C, 


labour ° is based on 
instead 
given 

two page discussion of nuclear reactors, a worked example 
and a discussion of 


dealing with a fluid catalytic cracker 


time lag in stirred vessels has been included. 


It is gratifying that the authors have not allowed this 
book to of date We 


edition of Volumes II and III with great interest. 


useful get out await the second 


J. F. 


DAVIDSON 


for hydrocarbons (19 pages), and for Table D (9 pages) 


in which a survey is given of representative cyclic hydro- 
carbons. 


In an appendix correction, figures for boiling points are 
given and nomographs for correcting observed boiling 
points at a certain pressure to-e.g 760 mm of mercury. 


The writer of this review expresses his gratitude for 
the enormous amount of work which Mr. S. W. Ferris 
must have undertaken to compile the data in question. 
The book will prove to very be useful. Analogous works 


are desirable in which other physical constants, e.g. 


viscosity and surface tension, can be considered, if possible 
at different temperatures. Perhaps the author can later 


find time to handle these constants also. 
H. I. 


WATERMAN 
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large scale produc if greet 


under controlle nditions 
on of farming 
therefore 


enthusiasn 


bie for large scalt 


ther of the 
Chlorella 


sorne mar 


superior 
grows in suspension U 


livision. It is autotroy that is 


f organic origin; it uses visible 


substance from carbon dioxide nitrogen 


compounds (such as nitrates or ammonia gether with 


phosphate, potassium and a number of other necessary 


ncluding the in pertect aracternized 


inorganic nutrients 


micro-nutrients or trace elements The dr ells contain 


sorne 0 protein and If fat, the ret r being mostly 


protein conta al t imino- 


“arbohvdrate The 


acids known to be essential in the ul t ithougt 


like all vegetable protein, it 
Feeding tests on rats 


omparable to yeast 


Primarily, chiorella is t 


provider of protein, which is 


in the world. The enthusiasts suggest that at any rate 


where 4yl tural iand is 


in some parts of the world 


scarce or over-populated - for instance 


Japan rT 


Israel —the scientifi 


justified by the high yi nit area of 


ground which is made possible by conti: is as against 


seasonal growth. Under such circumstances the process 


would not necessarily have to be economica competitive 
with conventional agriculture, as it would, for instance, in 
S.A. The sceptic might 
reply that only a virtually starving population would eat 
Chlorella, and that the effort might 
profitable be devoted to improvements in 


still 


order to be worth while in the I 


money and more 


agricultural 


practice, for which there is great scope in most 


countries. 


A plant for the large-scale culture of algae would be, in 


of the 
nutnment 


tank reactor 
The 


which ts drawn 


effect, a stirred type so 


{familiar to chemical engineers 


continuous 
medium is 
ff at the 


maintain 


culture 
itself must be ritat te 


and t 


fed continually to the 


rate. The culture 


same 


uniformity of composition keep the cells in suspen 


sion. The composition and feed-rate of the solution should 


be such as to maintain the optimum nutrient composition 


and cell concentration in the culture, which should be 


kept saturated with CQO, at an appropriate partial pressure 


(about 0-05 atm) and maintained at a suitable temperature 


heating or cooling may be required, depending on the 


limate It seems probable that satisfactory results can 


ynly be obtained by enclosing the culture under a trans- 


parent cover to keep in CO,., prevent evaporation (in 


regions where water is scarce) and exclude contaminating 


organisms 4 pilot plant constructed on these lines was 


constructed and operated from July to December 1951, in 


irculated through a wide 
of about 


Massachusetts. The culture was « 


polythene tube, which flattened out to a width 


4 ft and 


which was filled with 5 CO, in air 


contained an inch or so of liquid and a gas space 
Nearly 100 Ib (dry 
was and harvested, and the 


weight) of Chiorella grown 


results suggested that a vield of some 17$ tons per acre 
7:3 kg/m*) per year might be obtained from a plant of 
this type. This represents a yield of protein much higher 
than can be obtained from conventional crops, but it is 
clear that the capital charges and running costs (much 
of the 
a different 


land 


latter for power consumed in pumping) would be of 


rder from those for an equal area of agricu!tural 


There here to discuss the many technical 


arising in the design and operation of a con 


juous culture plant. These include the economic har 
the question of whether 
but the 


the danger of contamination 


vesting and drying of the algae 
the medium can be recycled with no treatment 
addition of fresh nutrients ; 
the supply of trace elements and the use 
buffer their 


unproductive absorption in the upper 


by protozoa ; 
f chelating agents to concentration ; the 
of light by 
and the possibility of mitigating it by making use 
effect 
The report contains many studies of such problems 
U.S.A., England 


waste 
layers 
of the 
others 


intermittent illumination and many 


by various contributors in the Japan, 


Israel and Germany 


It is difficult to come to any confident conclusion about 
the future of the large scale cultivation of algae. How- 
ever, two points are suggested by a study of this report. 
Firstly, in any project of this kind, competent chemical 
engineers should be brought in at an early stage if the 
of a 


matter is not to remain indefinitely at the level 


laboratory investigation; and secondly, the dangerous 
fascination of the technical problems involved may lead 
to the diversion of much valuable scientific and technical 


talent from work of more immediate promise. 


P. V. DANCKWERTS. 





